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ABSTRACT 

We  study  Weyl-Heisenberg  (=Gabor)  expansions  for  either  I/2(]R*^)  or  a  subspace  of  it.  These 
are  expansions  in  terms  of  the  spanning  set, 

X  =  {E^M^ip:keK,  leL.ipE^), 

where  K  and  L  are  some  discrete  lattices  in  IR^,  $  C  I/2(IR^)  is  finite,  E  is  the  translation 
operator,  and  M  is  a  modulation  operator.  Such  sets  X  are  known  as  WH  systems.  The  analysis 
of  the  ‘‘basis”  properties  of  WH  systems  (e.g.  being  a  frame  or  a  Riesz  basis)  is  our  central  topic, 
with  the  fiberization-decomposition  techniques  of  shift-invariant  systems,  developed  in  a  previous 
paper  of  us,  being  the  main  tool. 

Of  particular  interest  is  the  notion  of  the  adjoint  of  a  WH  set,  and  the  duality  principle  which 
characterizes  a  WH  (tight)  frame  in  term  of  the  stability  (orthonormality)  of  its  adjoint.  The 
actions  of  passing  to  the  adjoint  and  passing  to  the  dual  system  commute,  hence  the  dual  WH 
frame  can  be  computed  via  the  dual  basis  of  the  adjoint. 

Estimates  for  the  underlying  frame/basis  bounds  are  obtained  by  two  different  methods.  The 
Gramian  analysis  applies  to  all  WH  systems,  albeit  provides  estimates  that  might  be  quite  crude. 
This  approach  is  invoked  to  show  how,  under  only  mild  conditions  on  X,  a  frame  can  be  obtained 
by  oversampling  a  Bessel  sequence.  Finally,  finer  estimates  of  the  frame  bounds,  based  on  the  Zak 
transform,  are  obtained  for  a  large  collection  of  WH  systems. 
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Weyl-Heisenberg  frames  and  Riesz  bases  in  L2(]R‘^) 
Amos  Ron  and  Zuowei  Shen 


1.  Introduction 

1.1.  Frames,  Riesz  bases,  and  their  dual  systems 

The  present  paper  is  the  second  in  a  series  of  three,  all  devoted  to  the  study  of  shift-invariant 
frames  and  shift-invariant  stable  (=Riesz)  bases  for  H  :=  L2(1R.*^),  d  >  1,  or  a  subspace  of  it.  In 
the  first  paper,  [RSI],  we  studied  such  bases  under  the  mere  assumption  that  the  basis  set  can  be 
written  as  a  collection  of  shifts  (namely,  integer  translates)  of  a  set  of  generators  The  present 
paper  analyses  Weyl-Heisenberg  (=:WH,  known  also  as  Gaborian)  frames  and  stable  bases.  Aside 
from  specializing  the  general  methods  and  results  of  [RSI]  to  this  important  case,  we  exploit  here 
the  special  structure  of  the  WH  set,  and  in  particular  the  duality  between  the  shift  operator  and 
the  modulation  operator,  the  latter  being  absent  in  the  context  of  general  shift-invariant  sets.  In 
the  third  paper,  [RS3],  we  present  applications  of  the  results  of  [RSI]  to  wavelet  (or  affine)  frames. 
The  flavour  of  the  results  there  is  quite  different:  wavelet  sets  are  not  shift-invariant,  and  the  main 
effort  of  [RS3]  is  to  show  that,  nevertheless,  the  basic  analysis  of  [RSI]  does  apply  to  that  case  as 
well. 

Let  X  C  Z-2(1R‘^).  We  consider  A  as  a  possible  “basis”  set  for  L2(ffl‘^)i  or  for  some  closed 
subspace  of  it.  The  various  notions  of  a  “basis  set”  are  conveniently  defined  with  the  aid  of  the 
so-called  synthesis  operator  or  reconstruction  operator  T  :=  Tx  defined  by 

Tx  :  fo(^)  L2(1R'^)  ■■ 

Here,  £o(X)  is  the  collection  of  all  finitely  supported  sequences  in  ^2 (A).  If  T  is  bounded,  it  is 
extended  by  continuity  to  all  of  £2 (A).  We  use  the  notation  T  for  this  extension,  as  well. 

Definition  1.1.  A  is  a  basis  whenever  T  is  1-1  on  its  domain.  A  is  fundamental  if  ranT  is 
dense  in  L2(IR-‘^)-  If  T  is  bounded,  A  is  a  Bessel  set.  If  T  is  bounded  and  ranT  is  closed,  A  is  a 
frame.  Finally,  a  frame  which  is  also  a  basis  is  known  as  a  Riesz  (or  stable)  basis. 

Remark.  Some  of  the  articles  that  deal  with  frames  for  L2(iR'))  reserve  the  notion  of  “frame”  only 
to  the  case  that  we  refer  to  here  as  a  “fundamental  frame” . 

Note  that,  if  A  is  a  Riesz  basis,  T  has  a  bounded  inverse  which  acts  from  ranT  onto  ^2(A). 
We  denote  that  inverse  by  T~^.  If  A  is  not  a  Riesz  basis,  but  is  still  a  frame,  a  bounded  inverse 
still  exists  but  acts  from  ranT  onto  (kerT)-^.  This  pseudo-inverse  is  denoted  hereafter  by  T|“^, 
and  is  referred  to  as  the  partial  inverse  of  T. 

Another  way  to  define  the  above  “basis”  notions,  is  with  the  aid  of  the  analysis  operator  or 
the  decomposition  operator,  which  is  the  formal  adjoint  T*  :=  Tj^  of  T,  and  is  defined  by 

T*  :  L2{1R^)  ->  e2iX)  :  f  ^  {{f,x)Ux. 
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The  equivalent  definitions  via  the  adjoint  map  are  entirely  analogous:  the  definitions  of  frames  and 
Bessel  sets  remain  unchanged  (just  replace  T  by  T*).  The  fundament ality  of  X  amounts  to  the 
injectivity  of  T*,  and  a  Riesz  basis  is  a  frame  whose  corresponding  T*  is  surjective. 

Since  usually  the  injectivity  of  a  map  is  easier  to  check  than  its  surjectivity,  fundamental  frames 
are  usually  studied  via  T*,  while  Riesz  bases  axe  usually  analyzed  via  T.  That  tradition,  however, 
will  not  be  followed  in  the  present  paper,  since  our  techniques  are  critically  based  on  a  simultaneous 
analysis  of  both  the  decomposition  and  reconstruction  operators. 

In  order  to  decompose  and  reconstruct  functions  in  1/2 one  needs  to  complement  the 
given  frame  X  by  another  one,  the  so-called  dual  frame  which  is  defined  as 

RX, 

with 

The  map 

/  ^  TxT^xf  “  (/?  ^)  ^ 

is  then  an  orthogonal  projector  onto  ranTx?  hence  is  the  identity  on  L2(]R  )  in  case  the  frame 
X  is  fundamental.  In  this  regard,  computing  T^xf  decomposition  of  /,  and  computing 

TxT^^f  (from  the  given  sequence  f)  is  the  reconstruction  of  /.  The  operator  norms  HTxll 
and  ure  known  as  the  frame  bounds  of  X.  A  frame  whose  two  frame  bounds  coincide 

is  a  tight  frame,  and  it  is  well-known  that  the  dual  frame  of  a  tight  frame  is  (up  to  a  multiplicative 
constant)  the  original  frame  X  (i.e.,  R  =  cld  on  ranT). 

Frames  were  introduced  in  [DS]  in  the  context  of  non-harmonic  Fourier  series.  Frames  for 
L2  (IR)  were  quite  extensively  studied  in  the  literature,  with  the  focus  being  on  two  special  choices: 
wavelet  (or  affine)  frames^  and  Weyl-Heisenberg  (or  short-windowed  Fourier  transform)  frames.  We 
refer  to  the  surveys  [BW],  [HW]  and  [Dl],  the  books  [C]  and  [D2],  and  to  the  references  therein 
for  discussions  of  frames,  and  in  particular  wavelet  and  Weyl-Heisenberg  frames. 

1,2.  Weyl-Heisenberg  systems  defined 

A  Weyl-Heisenberg  (WH)  system  is  defined  here  with  respect  to  a  pair  of  lattices.  Here,  a 
lattice  (or,  more  precisely,  a  d-dimensional  lattice)  K  is  the  image  Ak'ZL  of  a  linear  invertible 
map  Ak  •  IR^  IR^-  The  column  vectors  of  Ak  generate  K.  The  volume  of  the  lattice 

\K\  ~  \det  Ak\ 

measures  the  sparsity  of  it,  and  depends  only  on  K.  The  dual  lattice  X  of  X  is  the  lattice  defined 
by 

X  :=  {/  €  IR^  :  /  •  fc  G  27rZ,  VA:  E  X}. 

Note  that  we  always  have 

(1.2)  \K\  \K\  =  {2nf. 
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Given  a  d-dimensional  lattice  K,  we  let 

stand  for  any  fundamental  domain  for  K,  i.e.,  a  set  for  which  Ufce/<-(fe  -1-  Qk)  is  an  (essential) 
partition  of  IR'^.  A  standard  choice  for  Qk  is  the  parallelepiped  Aa-[0  . .  1]*^. 

We  now  turn  to  the  definition  of  a  WH  system.  First,  for  t  €  H®*,  let  denote  the  translation 
operator  and  M*  the  modulation  operator,  i.e., 

+ 1),  M^:f^  etf,  teJR.\  fe  LjCIR'^), 

with  et  the  exponential  function 

Further,  let  $  C  I-2(IR‘^)  be  finite  (though  our  analysis  applies  to  infinite  $  as  well,  it  suffices  for 
all  practical  purposes  to  assume  that  $  is  finite).  Finally,  let  K,L  be  two  lattices  in  IR**.  Then,  we 
call  the  set 

X  =  {K,LU  ■■=  {E^M^<p  :  (p  e  k  E  K,  I  e  L} 

a  Weyl-Heisenberg  system  generated  by  $.  The  WH  system  is  normalized  whenever  ||<^|1  =  1, 
(pe^.  Also,  whenever  (p  is  a  singleton,  we  refer  to  A  as  a  principal  Weyl-Heisenberg  system 
(PWH  system,  for  short).  The  number 

(2^)^ 

is  the  density  parameter  den  {K,L)  of  {K,L),  and  is  also  referred  to  as  the  density 

denX 

of  a  PWH  system  X.  We  call  a  PWH  system  X  a  high-density  system  if  denX  >  1  and  a  low 
density  system  if  denX  <  1. 

1.3.  Layout  of  the  paper 

In  §2,  we  introduce  the  notion  of  the  adjoint  X*  of  a  PWH  system  X.  Our  duality  principle 
of  WH  sets  (Theorems  2.2  and  2.3)  deals  with  the  basic  relations  between  a  system,  its  adjoint, 
and  their  dual  systems.  In  general  terms,  the  duality  principle  exhibits  an  intimate  relation  be¬ 
tween  Tx>  and  (hence  also  between  Tx  and  TJ^,,  since  X**  =  X),  which  “almost”  says  that 
Tx>  =  cUiTxU2i  Ui,U2  unitary,  and  is  almost  as  useful  as  such  connection.  Indeed,  a  variety  of 
applications  of  the  duality  principle  are  then  presented;  some  generalize  known  univariate  observa¬ 
tions,  and  others  are  new  even  for  univariate  systems.  As  an  example  for  the  former,  a  fairly  trivial 
argument  is  invoked  to  show  (Corollary  2.10)  that  low-density  frames  are  never  fundamental.  We 
are  confident  that  further  applications  of  that  principle,  above  and  beyond  what  is  collected  in  the 
present  paper,  will  be  found  in  the  future,  and  that  the  duality  principle  will  become  one  of  the 
cornerstones  of  the  theory  of  WH  systems. 

In  the  first  part  of  §3  (§3. 1-3.3),  we  discuss  the  duality  principle  in  terms  of  the  fiberization 
techniques  of  [RSI].  In  §3.1,  our  basic  observation  concerning  that  principle  is  presented  (cf.  (3.1)), 
and  is  followed  by  the  proofs  of  Theorems  2.2  and  2.3.  Partial  unitary  relations  between  the 
operators  of  X  and  these  of  X*  are  discussed  in  §3.2. 
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In  §3.4,  we  derive  various  estimates  on  the  frame/Riesz  bounds  of  a  given  system.  These 
estimates  are  obtained  by  elementary  matrix-theory-manipulations  applied  to  the  Gramian  and 
dual  Gramian  matrices,  the  latter  are  the  basic  tools  of  the  shift-invariant  methods  of  [RSI]. 
The  estimates  are  related  to  the  estimates  of  [Dl]  (for  univariate  systems),  and  several  other 
references  (such  as  [T02]),  and  a  discussion  comparing  our  estimates  to  their  literature  counterparts 
is  included. 

Finally,  in  §4,  we  employ  the  Zak  transform  for  the  analysis  of  a  special  (still  large)  class  of  WH 
systems,  termed  in  this  paper  compressible  WH  systems,  and  which  are  reduced  in  the  univariate 
case  to  WH  systems  with  rational  density  parameter.  The  crux  here  is  the  observation  that  the 
infinite  Gramian  and  dual  Gramian  matrices  of  a  compressible  WH  sets  are  actually  matrix-valued 
convolution  operators,  with  the  relevant  matrix  being  of  finite  order.  The  Zak  transform  enters  the 
discussion  as  the  symbol  of  these  convolution  operators. 


1.4.  Acknowledgments  and  general  remarks 


Our  original  interest  in  shift-invariant  bases  stems  from  the  role  of  such  systems  in  Approx¬ 
imation  Theory  (e.g..  Box  Splines).  There,  non-fundamental  bases  are  the  rule  rather  than  the 
exception,  and  this  explains  our  genuine  interest  in  and  emphasis  on  non-fundamental  sets.  On  the 
other  hand,  frames  have  hardly  been  considered  in  Approximation  Theory  as  an  object  of  interest, 
and,  in  fact,  our  initial  development  of  the  frame  material  in  [RSl-3]  was  done  “from  scratch  . 
While  this  somewhat  cavalier  approach  might  have  had  its  own  advantages,  it  also,  inevitably, 
resulted  in  the  re-invention  of  known  and  even  classical  results  (the  Zak  transform  and  the  [DGM] 
painless  construction  of  WH  frames  were  among  our  early  “innovations”).  Communications  we  had 
in  late  1992  with  Chris  Heil  had  helped  us  in  drawing  connections  between  our  work  the  rich  frame 
literature. 

Our  first  presentation  of  the  duality  principle  (in  Oberwolfach,  Summer  1993)  had  led  to 
several  very  useful  discussions  with  Hans  Feichtinger  and  with  some  of  his  Vienna  group  people. 
In  particular,  numerical  experiments  conducted  by  Qiu  Sigang  had  helped  us  in  correcting  the 
constant  that  appears  in  Theorem  2.3. 

Selected  results  from  §2  (such  as  Theorems  2.2  and  2.3)  and  §4  are  announced  in  [RS2].  After 
receiving  a  copy  of  [RS2],  Meir  Zibulski  from  the  Technion,  Israel,  had  brought  to  our  attention  the 
articles  [T01,2]  and  [ZZ].  The  former  establishes  the  univariate  equivalent  of  part  (f)  of  Theorem 
2.2,  while  the  latter  derives  the  univariate  equivalents  of  part  (a,b)  of  Theorem  4.14,  as  well  as  of 


some  other  results  from  §4. 

About  the  time  we  were  essentially  done  with  the  present  endeavour,  Ingrid  Daubechies  had 
brought  to  our  attention  Janssen’s  paper  [J]  and  her  joint  paper  with  Landau  [DL].  Both  papers 
deal  with  the  same  phenomenon  that  we  describe  in  §2  here.  Specifically,  both  contain  statements 
equivalent  to  our  Theorem  2.3,  but  under  the  additional  assumptions  that  the  underlying  frame  is 
univariate  and  fundamental.  Both  also  contain  results  equivalent  to  the  univariate  case  of  parts 
(a,e)  of  our  Theorem  2.2,  including  the  connections  between  the  frame  bounds  asserted  in  that 
theorem.  It  is  probably  correct  to  consider  the  three  articles  [J],  [DL],  and  ours  as  simultaneous 


and  independent”,  and  it  is  worth  mentioning  that  the  techniques  employed  in  these  papers  are 
quite  different;  [J]  amplifies  the  approach  of  Tolimieri  and  Orr,  [DL]  invokes  what  they  call  “the 
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Wexler-Raz  identity”,  while  our  development  follows  the  fiberization  techniques  of  [RSI].  We 
decided  to  abstain  from  expanding  our  paper  in  directions  that  may  be  suggested  from  the  reading 
of  [J]  and  [DL],  with  the  following  single  exception:  in  §3.3,  we  show  how  the  Wexler-Raz  identity 
can  be  observed  by  using  our  decomposition-fiberization  techniques. 

We  would  like  to  extend  our  thanks  to  all  the  people  whose  contribution  is  detailed  above. 

A  final  remark:  our  results  are  always  derived  in  a  multivariate  setup,  and  deal  with  systems 
which  are  not  necessarily  fundamental  in  L2(]R‘^).  It  is  probably  true  that  no  significant  simplifica¬ 
tion  of  arguments  would  have  occurred,  had  we  chosen  to  restrict  attention  to  univariate  systems. 
In  contrast,  the  treatment  of  non-fundamental  systems  seems  to  be  harder  than  their  fundamental 
counterparts,  at  least  from  the  standpoint  of  the  tools  we  borrowed  from  [RSI]  for  either  case. 


2.  The  duality  principle  and  some  of  its  applications 

We  start  our  discussion  here  with  the  introduction  of  the  (new)  notion  of  the  adjoint  of  the 
PWH  system  {K,L)^. 

Given  (p  e  L2(IR‘^),  we  associate  each  X  =  {K,L)y,  with  another  PWH  system,  denoted  by 
X*,  referred  to  hereafter  as  the  adjoint  system  of  X,  and  defined  by 

X*  :=  (L,X)^. 

We  also  refer  to  {L,K)  as  the  adjoint  {K,L)*  of  (X,L).  Note  that  the  density  parameter  of  X*  is 
reciprocal  to  that  of  X: 

(2.1)  den(X)den(X*)  =  1. 

The  system  X  is  said  to  be  self-adjoint  if  X  =  X*  (i.e.,  if  K  =  L).  Note  that  all  self-adjoint 
systems  are  of  density  1  (as  follows  at  once  from  (2.1)),  but  not  vice-versa,  unless  d  =  1,  as  shown 
by  the  following  example. 

Example.  Assume  that  X  is  a  univariate  PWH  system.  Then,  (K,  L)  =  {pTL,  2'Kq7L)  for  some 
parameters  p, g  >  0.  Here,  den(X)  =  The  adjoint  system  is  {L,K)  =  {7Zfq,2T:7Zi/p), 

and  its  density  is,  indeed,  pq.  If  pq  =  1,  the  adjoint  system  coincides  with  the  original  system. 
Consequently,  a  univariate  PWH  system  is  self-adjoint  iff  its  density  is  1. 

We  prove  in  §3  the  following  result  concerning  the  connection  between  a  PWH  system  X  and 
its  adjoint  X*. 

Theorem  2.2.  Let  X  be  a  normalized  PWH  system.  Then: 

(a)  X  is  a  Bessel  system  if  and  only  if  X*  is  one.  In  that  case,  ||Tx*  |P  =  den(X*)||rxlp. 

(b)  Suppose  that  X  is  a  Bessel  system.  Then  X*  is  a  basis  if  and  only  if  X  is  fundamental. 

(c)  X*  is  a  frame  if  and  only  if  X  is  a  frame.  In  that  case,  =  den(X*)l|rx|“^ll”^. 

(d)  X*  is  a  tight  frame  if  and  only  if  X  is  a  tight  frame. 

(e)  X*  is  a  Riesz  basis  if  and  only  if  X  is  a  fundamental  frame. 

(f)  X*  is  an  orthonormal  basis  if  and  only  if  X  is  a  fundamental  tight  frame. 
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We  note  that  the  core  claims  in  the  above  theorem  are  (a-c),  with  the  rest  being  simple 
corollaries.  Specifically,  (d)  follows  directly  from  (a)+(c),  (e)  follows  from  (b)4-(c),  and  (f)  follows 
from  (d)+(e).  The  statement  in  (b)  is  actually  valid  even  when  X  is  not  Bessel,  if  one  defines  then 
appropriately  the  notion  of  a  “basis”.  Finally,  note  that  the  roles  of  X  and  X*  in  the  theorem  can 
certainly  be  interchanged  (since  X**  =  X). 

Remark.  The  paper  [T02]  precedes  us  in  observing  the  connection  (in  the  univariate  setup) 
between  a  set  (p2Z,27rgK)^  and  (K/g,27ra/p)^.  The  univariate  case  of  part  (f)  of  Theorem  2.2 
can  already  be  found  in  that  reference. 

The  actions  of  taking  adjoint  and  taking  dual  commute: 

Theorem  2.3.  Let  X  be  a  PWH  frame  generated  by  ip,  X*  its  adjoint  frame.  Then  the  dual  of 
the  adjoint  is  the  adjoint  of  its  dual;  more  precisely, 

=  (denX)(RxX)*. 


In  particular, 

=  (denX)RxV^- 

In  words,  the  generator  of  the  dual  frame  of  X* ,  is  the  same  as  the  generator  of  the  dual  frame  of 
X,  up  to  the  (important)  multiplicative  constant  denX 

It  is  easy  to  see,  and  is  well-known  (cf.  Theorem  3.2  of  [BW])  that  dual  frame  of  X  (K,  L),p 
is  of  the  form  (K,  L)^ ,  for  some  V'  €  La ,  hence  the  entire  difficulty  in  computing  the  dual  frame  of  X 
lies  in  computing  this  generator  This  is  particularly  diflficult,  since  no  bi-orthogonality  relations 
exist  between  X  and  its  dual  frame.  On  the  other  hand,  if  X  is  a  fundamental  frame,  its  adjoint 
is  a  Riesz  basis,  hence  the  dual  of  this  adjoint  is  characterized  by  the  standard  bi-orthogonahty 
relations.  Theorem  2.3  asserts  that  the  generator  of  the  dual  frame  of  X*  is,  up  to  the  multiplicative 
constant  denX,  the  same  as  the  generator  of  the  dual  frame  of  X,  and  thus  suggests  a  simpler 
avenue  for  computing  the  generator  of  the  dual  frame  of  the  fundamental  X. 

We  refer  to  the  relations  between  X  and  X*  that  are  expressed  in  Theorems  2.2  and  2.3  (as 
well  as  to  (3.1),  which  is  the  technical  ground  for  all  these  connections)  as  the  duality  principle 
of  WH  systems. 

Remark.  In  [RS2],  Theorem  2.3  is  announced  with  a  (-n  incorrect)  multiplicative  constant 
(denX)^/^.  The  fact  that  our  constant  was  flawed  was  revealed  in  numerical  experiments  conducted 
by  Hans  Feichtinger  and  Qui  Sigang  (from  Vienna),  and  we  are  grateful  to  them  for  pointing  out 
this  fact  to  us. 

The  rest  of  this  section  is  devoted  to  various  corollaries  and  applications  of  the  duality  principle. 
Corollary  2.4.  A  self-adjoint  X  is  a  fundamental  frame  if  and  only  if  it  is  a  Riesz  basis. 

This  property  is  immediate  from  (e)  of  Theorem  2.2.  The  “only  iP  implication  is  well-known 
(cf.  Theorem  3.7  of  [BW]  which  contains  a  proof  of  that  implication  for  cartesian  lattices). 

Another  interesting  corollary  is  the  following  (the  univariate  case  can  be  found  in  p.  81  of 
[D2]). 
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Corollary  2.5.  Let  X  be  a  normalized  fundamental  PWH  frame.  Then  the  frame  bounds  A,B 
of  X  satisfy 

A  <  denX  <  B. 

In  particular,  the  frame  bound(s)  of  a  normalized  fundamental  PWH  tight  frame  is  denX. 

Proof.  By  (e)  of  Theorem  2.2,  the  adjoint  X*  of  X  is  a  Riesz  basis.  Denoting  by  Aq,  Bq 
the  Riesz  bounds  of  X*,  it  is  clear  that  Aq  <  1  <  Bo  (since  Tx-  maps  a  sequence  of  norm  1,  viz. 
the  ^-sequence,  to  the  function  <p,  whose  norm  is  1  by  assumption).  An  application  of  Theorem 
2.2,  parts  (a,c),  then  yields  the  desired  results  on  the  frame  bounds  of  X.  A 

Another  application  is  a  new  proof  (and  in  several  dimensions)  for  the  “painless  construction 
of  tight  WH  frames”  that  was  done  in  [DGM].  For  that,  suppose  that  we  want  to  construct  a  fun¬ 
damental  frame  for  L2(R‘^)  of  the  form  (X,!)^,  (p  compactly  supported.  By  the  duality  principle, 
this  is  equivalent  to  constructing  a  Riesz  basis  (L,K),p.  Since  <p  is  compactly  supported,  we  can 
choose  L  sufficiently  sparse  to  guarantee  that  the  sets 

I  -f  supp  (p,  I  G  L 

are  disjoint.  In  such  case,  the  task  of  ensuring  that  (L,  X)^  is  a  Riesz  basis  is  reduced  to  ensuring 
that  Y  :=  (efcV^)^^^  is  a  Riesz  basis.  Furthermore,  if  Y  is  only  a  frame,  then  X*  is  a  frame,  too, 
with  the  same  frame  bounds,  hence  X  is  a  frame,  though  not  a  fundamental  one. 

Corollary  2.6.  Let  X  =  (X,  L),^  be  a  normalized  PWH  system  generated  by  a  compactly  sup¬ 
ported  function  <p,  and  assume  that  the  L-shifts  of(p  have  pairwise  disjoint  supports.  Then  (X,  L),^ 
is  a  frame  if  and  only  if  the  set  Y  :=  (efc<p)^£^  is  a  frame.  The  frame  bounds  of  X  are  denX  times 
the  frame  bounds  ofY.  Furthermore,  X  is  fundamental  iffY  is  a  Riesz  basis,  and  X  is  fundamental 
and  tight  iff  Y  is  orthonormal. 

As  is  well-known  (and  not  hard  to  prove),  Y  above  is  a  Riesz  basis  if  and  only  if  the  function 

:=  ( ^  \V>{-  +  A:)P)^/^ 

keK 

is  bounded  above  and  below  (i.e.,  away  from  0).  Also,  (by  Theorem  2.2.7  of  [RSI];  see  also  §7  of 
[BW])  y  is  a  frame  iff  Ipx  is  bounded  above  and  below  on  its  support  aY,  and  the  frame  bounds 
are 

(2-7)  \^\\\fK\\\^(„Yp  l-^l 

This  yields,  in  view  of  Corollary  2.6,  that  the  frame  bounds  of  X  are 
(2-8) 

The  computation  of  the  dual  frame  is  straightforward  here.  Indeed,  the  orthogonality  between 
the  L-shifts  of  <p  implies  that  the  function  tp  that  generates  the  dual  frame  of  Y  generates  also  the 
entire  dual  frame  of  X*.  The  simplification  is  then  due  to  the  fact  that  computing  the  dual  frame 
of  y  is  relatively  easy  (cf.  [RSI].  The  special  case  of  a  Riesz  basis  is  well-known).  The  generator 
rp  of  the  dual  frame  of  Y  is  the  function 
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with  the  understanding  that  0/0  :=  0.  This  function  generates  the  dual  frame  of  T,  hence  the  dual 
frame  of  X* .  In  view  of  Theorem  2.3,  we  obtain  that  the  frame  dual  to  X  is  generated  by 


rP  :=  \L\-^ 


JP_ 

* 


We  summarize  these  observations  as  follows: 

Corollary  2.9.  Let  X  =  {K,L)^  be  a  normalized  PWH  (fundamental)  system,  and  assume  that 
the  L-shifts  of  cp  have  pairwise  disjoint  supports.  Then: 

(a)  X  is  a  frame  if  and  only  if  the  function  pK  is  bounded  above  and  away  from  0  on  its  support 
gY  (with  aY  =  IR.‘*  in  the  fundamental  case).  Furthermore,  the  frame  bounds  of  X  are 

l^lll^K|lL(c.y),  and  |£|||l/<?xllZ^(^y)- 

(b)  Assuming  X  to  be  sl  frame,  the  dual  frame  of  X  is  the  system  where 

yk 

The  analysis  is  applicable  also  in  case  (p  rather  than  ip  is  compactly  supported.  One  simply 
applies  the  same  arguments  to  the  system  [L^K)^  instead  of  (frT,  L)<^.  Thus,  for  this  case,  the 

X-shifts  of  (p  are  required  to  have  disjoint  supports. 

Another  interesting  result  which  follows  directly  from  the  duality  principle  and  Theorem  2.3 
is  the  following: 

Corollary  2.10,  Let  X  be  a  PWH  system.  Then,  X  can  never  be  a  fundamental  frame  for  L20fV^) 
unless  den  AT  >  1. 

The  (strikingly  simple)  proof  of  that  fact  is  based  on  the  following 
Lemma  2.11.  Let  ip  be  the  generator  of  a  PWH  fundamental  frame  X,  Let  ip  be  the  generator 
of  the  frame  dual  to  X,  i,e,  xp  =  Then 

=  (denX)"^ 

Proof.  By  (e)  of  Theorem  2.2,  X*  is  a  Riesz  basis,  whose  dual  system  is  generated  (Theo¬ 
rem  2.3)  by  (denX)V>.  The  bi-orthogonality  relations  between  a  Riesz  basis  and  its  dual  basis  then 
imply  that 

{ip,  {denX)ip)  =  1.  4 


Proof  of  Corollary  2.10.  Let  RX  be  the  frame  dual  to  X.  We  recall  that,  for  every  /  €  L2, 
is  the  (unique)  minimal-norm  sequence  in  the  pre-image  T^f.  Choosing  /  :=  p,  and  taking 
into  account  the  fact  that  contains  one  sequence  of  norm  1  (viz.  the  sequence  that  is  1  at 

p  and  zero  elsewhere),  we  conclude  that  ||T]axV^II  ^  L  particular,  with  ip  the  generator  of  the 
dual  frame,  \{p,ip)\  <  1  (since  this  number  is  one  of  the  values  of  the  sequence  Combining 

the  above  with  Lemma  2.11,  we  conclude  that 

(denX)-^  =  {p,ip)<l. 

Hence,  indeed,  X  is  a  high-density  system.  ^ 
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It  is  known  that  the  univariate  analog  of  Corollary  2.10  does  not  require  X  to  be  a  frame,  only 
a  fundamental  set.  However,  the  simple  argument  for  the  univariate  case  (that  appears  in  p.  978  of 
[Dl])  applies  only  to  the  case  when  denX  is  rational  (which  is  the  univariate  counterpart  of  what 
we  call  here  compressible  PWH  systems,  see  §4,  particularly  Corollary  4.16),  and  with  a  different 
argument  (that  invokes  Rieffel’s  results  from  [R])  required  for  a  general  PWH  system. 

The  argument  that  was  used  to  prove  Corollary  2.10  easily  implies  the  following  further  esti¬ 
mate. 

Proposition  2.12.  Let  Y  be  the  frame  dual  to  the  PWH  fundamental  frame  X.  Let  be  the 
generator  of  the  frame  X,  and  ip  the  generator  ofY.  Then, 

|(<^,y)P<l-(denX)-2. 


We  had  mentioned  before  the  basic  fact  that  a  frame  and  its  dual  system  fail  to  satisfy  the 
bi-orthogonality  relations  that  are  satisfied  in  the  Riesz  basis  case.  However,  for  certain  PWH 
systems  X,  the  bi-orthogonality  relations  are  preserved,  though  only  to  a  limited  extent.  This  is 
the  case  of  a  sup-adjoint  PWH  system,  defined  as  follows. 

Definition  2.13.  A  PWH  system  X  =  {K,L)^  is  called  sup-adjoint  if  it  contains  its  adjoint  X* . 
Equivalently,  \i  L  C.  K. 

Note  that  the  assumption  L  d  K  implies  that  K  C  L,  implying  thereby  that  X*  =  {L,  K)^  C 
{K,L)ip.  Thus,  indeed,  the  condition  X*  C  X  is  equivalent  to  Lc  K.  Note  that  denX  here  is  the 
index  of  L  in  K,  hence  is  an  integral  number. 

Example.  Assume  K  =  7L^.  Then  the  condition  Z  C  S'*  is  equivalent  to  the  condition  that  L  is 
a  superset  of  2tt7L'^.  Therefore,  (S'*,  L)^  is  sup-adjoint  if  and  only  if  L  is  superlattice  of  27rS'*.  In 
particular,  all  PWH  systems  of  the  form  (S'*,27rS'*/n),p,  n  integer,  are  sup-adjoint. 

Theorem  2.14.  Every  fundamental  sup-adjoint  (tight)  frame  X  :=  (K,  L)^  is  a  union  of  the  form 

iK,L)^=Ug^G{K,L); 


of  (denX)^  Riesz  (orthogonal)  bases.  Here,  G  C  {K,L)^. 

Proof.  Let  (Xjjjg/  be  the  cosets  of  X*  in  X,  and,  with  Y  the  dual  frame  of  X,  let  (Fi)* 
be  the  corresponding  cosets  of  Y*  in  Y.  By  Theorem  2.3,  (denX)y*  is  the  dual  frame  of  X*. 
Assuming  X  to  be  fundamental,  we  know  (from  (e)  of  Theorem  2.2)  that  the  pair  (X*,  (denX)y*) 
is  bi-orthogonal,  and  from  that  it  easily  follows  that  each  pair  (Xj,  (denX)yi)  also  consists  of  bi- 
orthogonal  Riesz  bases.  Since  each  Xj  is  of  the  form  g  +  X*,  g  &  X,  the  result  for  a  fundamental 
frame  follows,  while  the  result  for  a  fundamental  tight  frame  follows  from  the  orthonormality  of 
X*  (Theorem  2.2,  part  (f)).  ♦ 
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Example.  Suppose  that  X  =  {7Z‘^,2Tv7Z^/n)^,  with  n  positive  integer.  Then  X*  =  {n7L'^,2n7L'^)^, 
hence  X  is  indeed  sup-adjoint.  The  cosets  of  X*  in  X  are 

(fc  +  ttZZ^, //n  +  27rZ^)<^,  = 

A  necessary  and  sufficient  condition  for  X  to  be  a  fundamental  tight  frame  is  the  orthogonality 
of  the  subset  (nK^,27rZ'^)^.  In  particular,  we  cannot  get  a  sup-adjoint  fundamental  tight  frame 
whose  elements  are  not  partially  orthogonal  in  the  above  sense! 


While  fundamental  frames  and  fundamental  tight  frames  can  be  constructed  with  some  ease  if 
the  generator  (p  is  known  to  be  either  compactly  supported  or  band-limited,  such  frames  can  also 
be  generated  by  functions  which  are  neither  compactly  supported,  nor  band-limited.  However,  if 
the  generator  cp  or  its  Fourier  transform  (p  is  positive  (as  is  the  case  with  the  Gaussian  kernel)  one 
cannot  use  p  to  construct  a  fundamental  tight  frame.  In  fact,  a  slightly  more  general  result  is  true. 

Theorem  2.15.  Let  p  be  a  generator  of  a  fundamental  frame  X.  Let  'll)  be  the  generator  of  the 
dual  frame.  Then: 

(a)  If  p  >  0  a.e.,  then  ^  assumes  positive  and  negative  values  in  a  non-trivia]  way  (i.e.,  on  sets  of 
positive  measures^. 

(b)  If  p  >  0  a.e.,  then  'ip  assumes  positive  and  negative  values  in  a  non-trivial  way. 

Proof.  If  X  =  {K,  L)^  is  a  fundamental  frame,  then,  by  Theorem  2.2,  X*  is  a  Riesz  basis, 
and,  by  Theorem  2.3,  ip  is  perpendicular  to  X*\p.  Since  this  latter  set  contains  functions  that  are 
positive  a.e.  (of  the  form  0p,  I  G  IR^),  it  follows  that  ip  cannot  be  essentially  of  one  sign.  This 
proves  (a),  while  (b)  follows  by  switching  to  the  fundamental  frame  {L,K)^  and  invoking  (a).  A 

Since  the  dual  frame  of  a  fundamental  tight  frame  {K,  is  generated  by  a  scalar  multiple  of 
p,  we  conclude  the  following  from  the  above  theorem. 

Corollary  2.16.  Assume  p  G  L2,  and  either  p  >  0  a.e.  or  p  >  0  a.e.  Then,  there  exists  no  PWH 
fundamental  tight  frame  generated  by  p. 

As  Corollary  2.6  already  indicates,  the  last  result  cannot  be  extended  to  functions  p  which  are 
only  positive  on  their  support.  However,  the  argument  used  above  leads  to  the  following  partial 
converse  of  Corollary  2.6. 

Corollary  2.17.  Let  X  =  {K,  L)^  be  a  fundamental  tight  frame,  and  assume  that  p>0  frespec- 
tively,  p>  0)  a.e.  on  its  support.  Then  the  sets  {I  +  suppp)^^'^  (respectively,  (k-i’Snppp)^^^)  are 
essentially  disjoint. 

Proof.  Since  X  is  fundamental  and  tight,  its  adjoint  X*  is  orthogonal.  Thus,  if  suppcp^ 
(l-hsuppp)  has  a  positive  measure,  the  positivity  of  p  implies  that  (p,£l  V)  ^  hence  I  ^  L. 

The  argument  for  p  is  essentially  the  same.  ^ 
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Example.  If  supp^  =  [0 ..  a]  and  y?  >  0  on  (0,a),  then  Corollary  2.6  combined  with  Corollary 
2.17  show  that  {p7Z,q7Z)^  is  a  fundamental  tight  frame  if  and  only  if  the  following  two  conditions 
hold: 

(a)  aq  <  2it. 

(b)  T,je7z.  l¥’(-  +  iP)P  is  constant. 

Finally,  we  have  already  commented  on  the  value  of  Theorem  2.3  for  the  computation  of  the 
generator  of  the  dual  frame  of  a  fundamental  frame.  We  add  below  some  further  details  in  that 
direction. 

Theorem  2.3,  together  with  the  duality  principle,  allows  us  to  derive  the  following  characteri¬ 
zation  of  the  generator  r/)  of  the  dual  of  a  fundamental  PWH  frame. 

Corollary  2.18.  Let  X  =  iK,L)^  be  a  fundamental  frame.  Then  the  generator  ip  of  the  dual 
frame  of  X  is  the  only  function  in  L2(JR"‘)  that  satisfies  the  following  two  conditions: 

(a)  Ip  €  ranTx*. 

(b)  ZkeK  E'^ipW^  =  I  e  1. 

Proof.  Since  X  is  a  fundamental  frame,  its  adjoint  X*  is  a  Riesz  basis.  By  Theorem  2.3, 
(denX)^  generates  the  dual  basis  of  X*  =  (L,K)^.  Thus,  ip  is  characterized  by  the  condition 
Ip  €  ranTx*  together  with  the  bi-orthogonality  conditions 

{iP,  =  (denX)-i<y,,o4,o,  0,  k)  E  Lx  K. 

This  last  condition  is  equivalent  to  the  stated  condition  (b)  (after  talcing  into  account  the  fact  that 
(denX)-i  =  \L\-^\K\.)  4 

The  result  reduces  the  computation  of  ip  to  solving  a  bi-infinite  linear  system  of  equations. 
Under  more  restrictive  conditions  on  the  pair  {K,  L),  we  will  show  in  section  4  that  the  computation 
of  Ip  can  be  achieved  by  solving  a  finite  system  of  linear  equations  (with  function- valued  coefficients). 

For  a  fundamental  tight  frame,  the  last  corollary  (or,  more  directly,  (f)  of  Theorem  2.2)  is 
specialized  as  follows. 

Corollary  2.19.  Let  X  be  a  normalized  PWH  system  generated  by  <p.  Then  X  is  a  tight  frame 
for  L2(IR-‘*)  If  3.nd  only  if  the  following  conditions  hold: 

(2.20) 

keK 

The  case  discussed  in  Corollary  2.6  can  now  be  identified  as  a  simple  instance  where  condition 
(2.20)  automatically  holds  for  1^0. 
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3.  Gramian  Analysis 


3.1.  The  Gramian  matrices  and  the  duality  principle 

Let  X  =  (A,L)$  be  a  WH  system  ($,  say,  finite).  Then,  X  is  A-shift-invariant,  in  the  sense 
that  each  of  the  shift  operators 

E’^  f{-  +  k),  keK 

maps  X  M  onto  itself.  The  reference  [RSI]  suggests,  for  a  general  K-shift-invariant  X,  the  study  of 
Tx  and  via  a  decomposition  process  of  these  two  operators  into  a  collection  of  “fiber  operators 
that  are  simpler  for  analysis.  For  that  purpose,  we  consider  the  set  X  as  the  collection  of  all 
A-shifts  of  the  set 

:=  {ei<p  :  v?  €  $,  /  €  L}. 

The  basic  object  in  the  approach  of  [RSI]  was  the  pre- Gramian  matrix  J  :=  Jx-  In  the 
present  case  (i.e.,  of  the  WH  system  X),  the  pre-Gramian  is  an  infinite  matrix  with  L2-valued 
entries,  whose  rows  are  indexed  by  K,  whose  columns  are  indexed  by  and  whose  (A:,  G 

{K  X  L$)-entry  is 

\K\-^^^(pi-  +  k  +  l). 

We  let  J{w)  be  the  evaluation  of  J  at  any  w  6  IR*^.  ^ 

The  formal  adjoint  J*  of  J  is  the  matrix  indexed  by  x  K,  with  entries 

lX|-i/2^(.  + A:  +  /). 

We  use  the  collection  (J(w))u,£iRd  for  the  representation  of  Tx,  and  the  collection  (J*(n;)),„e]R<<  for 
the  representation  of  T^- 

Alrecidy  at  this  initial  stage,  we  are  ready  to  present  the  most  crucial  observation  concerning 
the  duality  principle.  With  X  the  PWH  set  (A,  L)^,  let  us  compute  the  «;-value  of  the  pre-Gramian 
Jx^iw)  of  the  adjoint  X*  =  (L,K)^:  the  entries  of  the  pre-Gramian  of  this  adjoint  are 

\L\-^/^ip{w  +  l  +  k),  il,k)€  (L,  K),  w  e 

It  follows,  thus,  that,  for  every  w  €  IR'*, 

(3.1)  Jx- H  =  =  (denXr'^^^J^iw). 

\L\y^ 

Thus,  roughly  speaking  (that  is,  ignoring  the  unitary  transformation  of  taking  complex  conjugation, 
and  ignoring  the  multiplicative  constant  (denX)"^/^),  the  pre-Gramian  of  X*  equals  pointwise  to 
the  adjoint  pre-Gramian  of  X;  this  is  the  essence  of  the  duality  principle. 

At  this  point,  the  discussion  can  be  advanced  in  two  different  complementary  ways.  The  first, 
which  we  present  here,  fully  invokes  the  fiberization  results  of  [RSlj.  Roughly  speaking,  these 
results  imply  that  the  basic  norm  bounds  of  Tx  can  be  computed  via  a  separate  inspection  of  each 
fiber  Jx{w)  of  Jx-  This,  when  combined  with  (3.1),  will  provide  an  immediate  proof  for  almost 
all  statements  in  Theorem  2.2.  The  other  approach  aims  at  direct  connections  between  the  two 
operators  Tx-  and  Tf^-,  and  is  presented  in  the  next  subsection. 
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Though  the  entire  study  of  the  decomposition  idea  can  be  performed  by  the  pre-Gramian 
and  its  adjoint,  we  found  it  more  convenient  to  express  the  various  results  in  terms  of  Hermitian 
matrices.  That  is,  the  matrix 

G{w)  :=  Gx{w)  :=  J*{w)J{w)  =  +  ^  +  +  ^  + 

k^K 


and 


G{w)\=Gx{w):=  J{w)J*{w)  =  \K\  (p{w-\-k  +  l)(p{w  +  k' + 


We  refer  to  G  as  the  Gramian  of  X  and  to  G  as  the  dual  Gramian  of  X.  The  entries  of 
both  matrices  can  be  easily  shown  to  be  locally  integrable,  hence,  in  particular,  are  well-defined 
a.e.  We  consider  each  (more  precisely,  almost  each)  G{u))  as  a  densely  defined  operator  from 
£2{L^)  into  itself,  and,  similarly,  each  G(u>)  as  a  densely  defined  operator  from  t^iK)  into  itself. 
Whenever  the  relevant  operator  is  well-defined  and  bounded,  it  is  extended  by  continuity  to  the 
entire  £2 (■£'$)  ihiK)).  The  inverse  operators  G(^^;)“^  G(^(;)“^  and  the  pseudo- inverse  operators 
G(to)j“^,  G(i/;)p^  are  defined  in  the  same  way  described  in  the  introduction. 

The  pertinent  result  here,  which  is  stated  below,  follows  firom  Corollary  3.2.2,  Theorem  3.3.5 
and  Theorem  3.4.1  of  [RSI]. 


Theorem  3.2.  Let  X  be  a  WH system  as  above,  associated  with  a  Gramian  G  and  a  dual  Gramian 
G.  Consider  the  following  functions  (if  G(w)  is  the  0~operator,  we  define  its  partial  inverse  to  be 
0,  too.  Also,  if  the  underlying  operator  is  not  well-defined  or  is  unbounded,  its  norm  equals  00,  by 
definition): 

Q 

G* 

Q- 
G*- 

^r 


Gx 

IR-j. 

:  w 

H4 

\\G{w)\\ 

Qx 

IR+ 

:  w 

II^HII 

IR*' 

]R-i- 

:  w 

1-4 

07 

IR'* 

:  w 

h4 

||G(«;)-i|| 

^\x 

IR-f 

:  w 

H4 

WGHi^W 

B'^ 

-> 

]R-f. 

:  w 

»-4 

l|G(u;)f'|| 

Then,  the  following  is  true. 

(a)  The  following  conditions  are  equivalent: 

(i)  X  is  a  Bessel  system. 

(ii)  G  €  Loo- 
(Hi)  G*  G  Loo- 

Furthermore,  ||rp  =  ||T-f  =  ||ff||i.  =  ||5-||i,.. 

(b)  Assume  X  is  a  Bessel  system.  Then  the  following  conditions  are  equivalent. 

(i)  X  is  a  frame. 

(ii)  G^  €  Loo- 

(iii)  G(~  e  Loo- 

Furthermore,  ||Tj-‘f  =  =  ||5|-||t.  =  ||e|-|U.. 

(c)  Assume  X  is  a  Bessel  system.  Then  the  following  conditions  are  equivalent: 
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(i)  X  is  a  Riesz  basis. 

(ii)  Q~  £  Loo- 

Furthermore,  ||T“^||^  =  ||^~||loo- 

(d)  Assume  X  is  a  Bessel  system.  Then  the  following  conditions  are  equivalent: 

(i)  X  is  a  fundamental  frame. 

(ii)  Q*-  e  Loo. 

Furthermore,  ||r*"^|P  = 

(e)  Assume  X  is  a  frame.  Then  the  following  conditions  are  equivalent: 

(i)  X  is  a  tight  frame. 

fii;a*=apa.e. 

(Hi)  Q  =  g^,  a.e. 

If  X  is  fundamental,  the  last  three  conditions  are  also  equivalent  to: 

(iv)  g*  =  g*-  a.e. 

(f)  Assume  X  is  a  Riesz  basis.  Then  the  following  conditions  are  equivalent. 

(i)  X  is  orthogonal. 

(ii)  G  =  Q~  a.e. 

Since  taking  complex  conjugation  is  certainly  a  unitary  operation,  we  conclude  from  (3.1)  that 

g^  =  {denX)gx»', 

g*x  ~  (den , 

etc.  Parts  (a,c,e)  of  Theorem  2.2  then  follow  from  these  relations  when  combined  with  the  various 
conditions  of  Theorem  3.2.  As  mentioned  before,  parts  (d,f)  of  Theorem  2.2  follow  from  (a,c,e), 
and,  finally,  a  direct  proof  of  (b)  is  given  in  §3.2. 

We  now  prove  Theorem  2.3.  The  basic  idea  behind  the  proof  is  relatively  simple:  if  Y  = 
(K,  L),p  is  the  dual  frame  of  X  =  (K,  L)^,  then  the  two  operators  T^Ty  and  TyT^  are  orthogonal 
projectors,  implying  that  the  corresponding  pre-Gramian  products  Jx^y  orthogonal 

projectors,  too.  At  this  point  our  decomposition  techniques  are  employed:  the  analysis  part  of 
Lemma  4.1  of  [RSI]  asserts  that  almost  each  of  the  corresponding  fiber  operators  is  an  orthogonal 
projector,  (3.1)  then  converts  that  property  to  the  fibers  of  the  corresponding  adjoint  systems,  and 
the  synthesis  part  of  that  Lemma  4.1  is  then  invoked  to  conclude  that  (up  to  the  right  multiplicative 
constant)  and  Ty.T^.  are  orthogonal  projectors.  However,  that,  in  general,  does  not  imply 

that  y *  is  the  dual  frame  of  X* ,  until  one  verifies  that  X*  and  Y*  span  the  same  subspace  of  L2(IR‘^) 
(cf.  Corollary  1.3.9  of  [RSI]).  The  fact  that  X*  and  Y*,  indeed,  span  the  same  space  follows  from 
the  fact  that  the  kernels  of  the  synthesis  operators  of  two  dual  systems  (viz.  X  and  Y)  are  identical 
(Proposition  1.3.7  of  [RSI]). 

However,  the  entire  argument  sketched  above,  including  all  missing  details,  had  already  been 
detailed  in  Corollary  4.4  of  [RSI],  which,  as  a  matter  of  fact,  was  tailored  there  for  the  proof  of 
Theorem  2.3.  All  we  need  here  is  to  carefully  compute  the  constants  that  arise  when  converting 
the  shift-invariant  setup  of  that  corollary  from  the  lattice  to  general  lattices. 
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Proof  of  Theorem  2,3.  With  =  Rx<p  the  generator  of  the  frame  dual  to  X,  set  :=  denX^, 
and  Y  :=  Our  objective  is  to  show  that  Y  is  the  dual  frame  of  X*, 

The  basis  of  the  duality  principle  is  the  identification  of  Jx*  3^  (denX)~’^^^  With  KX  the 
dual  frame  of  X,  our  definition  of  'ipo  and  Y  imply  also  that 

Jy  =  {denX)^^^J^x* 

This  establishes  the  relation  we  need  for  the  application  of  Corollary  4.4  of  [RSI],  with  the  following 
additional  remark:  the  corollary  assumes  the  underlying  systems  to  be  Z'^-shift-invariant,  hence 
requires  relations  such  as  ^x^kx  —  Jx-^Y-  Though  our  systems  are  K-  or  L-shift  invariant,  the 
same  relation  should  still  be  required  (i.e.,  without  any  normalization  factors;  this  is  due  to  the  fact 
that  we  use  throughout  normalized  pre-Gramians).  We  had  normalized  V’o  exactly  for  that  purpose, 
and  consequently  obtain  the  relations  “just  right”  (save  an  extra  harmless  complex-conjugation). 

With  these  clarification  in  hand,  we  invoke  the  above-mentioned  Corollary  4.4  to  obtain  that, 
indeed,  Y  is  the  dual  frame  of  X*.  4 

3.2.  Unitary  relations  between  the  analysis  operator  of  a  system  and  the  synthesis 
operator  of  its  adjoint 

We  have  not  discussed  in  the  previous  subsection  the  exact  meaning  of  the  presentation  pro¬ 
vided  by  (J(to))w  to  T  (and,  analogously,  by  {J*{w))w  to  T*).  A  complete  discussion  of  that  point 
(in  a  more  general  setup),  including  the  issue  of  well-definedness  of  these  operators  can  be  found  in 
§1  of  [RSI].  So  far,  we  had  circumvented  entirely  that  point,  and  invoked  instead  the  ready-to-use 
Theorem  3.2.  However,  for  the  sake  of  the  present  development,  we  need  at  least  a  basic  grasp 
of  that  representation.  All  WH  systems  considered  in  the  present  subsection  are  assumed  to  be 
Bessel. 

Given  an  indexed  set  if,  and  an  open  set  Q.  C  11“^,  let 

L2{^,K) 

be  the  Hilbert  space  of  all  functions  from  9.  to  t^iK)  which  are  measurable  and  square-integrable 
(in  the  sense  of  [H]).  Given  a  WH  system  X  =  we  may  identify  the  space  4  (A)  with 

1^2(0^,  L)  with  the  aid  of  the  unitary  map 

£2{X)  9  c  i->  U\C, 


where 

(3.3)  Uic{w,  1)  :=  (w). 

Here,  is  the  Fourier  series  of  the  restriction  of  c  to  {K,l),  and  is  any  fundamental  domain 
for  K  (i.e.,  9~  -I-  AT  is  an  essential  partition  of  IR**).  We  also  identify  L2(IR‘^)  with  L2{9~,  K)  via 
the  map 

(3.4)  L2(]R‘^)  9  /  f/2/, 
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with  U2f{w,k)  :=  f{w  +  k).  Interpreting  Jx  as  the  operator 

Jx  :  L2{^^^L)  — >  L2{^j^yK)^ 

defined  by  {Jx'^){'^r)  *=  it  follows,  [RSI],  that 

(3.5)  U2TX  =  JxUi, 
hence  that  also 

(3.6)  UiT*^  =  rxU2^ 

However,  it  will  be  erroneous  to  conclude,  in  view  of  the  above  and  (3.1),  that  we  had  established 
an  operator  relation  between,  say,  Tx*  and  T^.  The  reason  is  that,  while,  up  to  a  constant,  (3.1) 
indeed  shows  that  both  and  Jx*  are  synthesized  from  the  same  fiber  matrices,  the  two  operators 
differ  in  their  domain  and  target  spaces:  while  acts  from  L2{^j^^K)  into  X2(^^^?^)?  Jx*  acts 
from  L2{flL,K)  into  12(^1,  L).  Thus,  finding  exact  connections  between  the  relevant  operators  is 
more  subtle  than  it  may  look  like. 

Before  we  proceed,  we  sidetrack  to  list  one  corollary  of  the  above  discussion  of  independent 
interest  (though  quite  unrelated  to  the  present  course  of  development):  in  the  form  described  above, 
it  is  clear  that  every  space  L2{^^K)  C  L2{^j^^K)  is  an  invariant  subspace  of  the  operator  JyJx'i 
for  any  Bessel  X  =  {K,L)^  and  Y  =  (if,L)^.  Since  the  space  L2{^,K)  represents  all  functions 
whose  Fourier  transform  is  supported  on  +  iiT,  we  arrive  at  the  following  conclusion. 

Corollary  3.7.  Let  X  =  {K,L)^  and  Y  =  {K,L)^  be  any  WH  Bessel  systems.  Let  a  be  a 
measurable  subset  of  IR^.  Then: 

(a)  If  a  is  invariant  under  K -shifts,  then  the  space  of  all  L2’’functions  whose  Fourier  transform  is 
supported  on  a  is  an  invariant  subspace  ofTyT^- 

(b)  If  a  is  invariant  under  L-shifts,  then  the  space  of  all  L2-functions  supported  on  a  is  an  invariant 
subspace  of  TyT^ . 

Proof,  The  first  claim  was  proved  in  the  paragraph  preceding  the  corollary  (the  K- 
invariance  of  a  implies  the  existence  of  C  with  0  =  +  K).  The  second  claim  follows 

by  an  application  of  the  first  to  X  :=  (L,  K)^,  ^ 

It  will  be  technically  more  convenient  to  state  our  subsequent  results  in  terms  of  the  modified 
adjoint  _ 

X'  :=  (Z,ff)^(-..)- 

The  reason  is  that  Jx'  =  Jx*^  hence  the  connection  between  Jx  Jx*  is  somewhat  nicer 
than  (3.1): 

(3.8) 

This  implies  that  the  details  of  unitary  connections  that  we  are  going  to  discuss  below  are  slightly 
simplified  if  we  replace  X*  by  X'. 

We  have  just  mentioned  the  main  difficulty  in  establishing  operator  identities  that  connect  the 
X  and  X*  (or  X').  This  difficulty  disappears  if  we  assume  X  to  be  self-adjoint.  Indeed,  in  this 
case  L  =  K,  hence  the  pre-Gramian  identities  (3.1),  (3.8)  can  be  interpreted  as  identities  between 
operators.  Thus,  in  the  self-adjoint  case,  (3.1),  (3.5)  and  (3.6)  imply  the  following: 
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Theorem  3.9.  Let  X  be  a  self-adjoint  WH  system.  Then  there  exists  an  isometry  U  :  L2(IR‘^) 
such  that 

T;c  =  UTx'U. 

Consequently,  the  operators  T^Tx  and  are  unitarily  equivalent. 

In  fact,  the  previous  discussion  reveals  a  possible  operator  U.  It  is  composed  of  the  Fourier 
transform,  followed  by  TP,  when  y  =  {K,L)x,  and  x  the  characteristic  function  of  = 

In  more  genercil  setups,  it  is  easier  to  derive  connections  between  Tp  and  Tx'  if  one  is  will¬ 
ing  to  restrict  the  domain  of  these  operators.  Specifically,  let  fi  be  any  measurable  subset  of 
Q~  n  Treating  L2(fi,  C)  (with  C  either  K  or  L)  as  a  subspace  of  ^),  and  recalling  the 

isometry  between  this  latter  space  and  L2{1R^)  (cf.  (3.4)),  one  can  find  a  subspace  H  :=  Hax 
of  L2(IR‘^)  (which  is  the  space  of  all  functions  whose  Fourier  transform  is  supported  on  +  £) 
such  that  U2,cH  =  L2{fl,C),  with  U2,c  defined  similarly  to  (3.4).  In  the  same  manner,  following 
on  the  definition  of  Ui,  (3.3),  we  may  find  subspaces  of  t2{X')  and  ^  of  f2{X),  together 
with  associated  unitary  maps  C/i,/,  (respectively  ~)  that  map  (respectively,  S^  ~)  onto 

L2{fl,K)  C  L2(flL,K)  (respectively,  L2{fl,L)  C  L2{fl^,L)).  Finally,  we  denote  by  I~  the  “iden¬ 
tity  map”  that  transforms  L2{fl,K)  from  a  subspace  of  L2{fl'^,K)  to  a  subspace  of  L2{flL,K). 
Thus  we  have  the  chain  of  isometries 


H. 


il,K  u  J 


L2{fl,K) 


2,K 


L2{fl,K) 


Sn,L- 


Setting 


we  conclude  that  V~  is  an  isometry  between  ^  and  In  a  similar  fashion,  one  can  construct 
an  isometry  Vl  =  ^  between  and 

The  key  point  in  all  this  development  is  that,  on  L2{fl,K),  (3.1)  gives  rise  to  the  rigorous 
operator  relation 

(denX)  ^  =  IiJx'^K' 

Therefore,  we  may  combine  the  relations  (3.5),  (3.6)  (applied  to  X  and  X')  to  conclude  that  on 


(denX)-i/2r*  =  (denX)-i/2c/* ~ 


U*  ~I  J  I~U  ~  =  v  U*  J  U  v~  =  V  T 


K' 


We  have,  thus,  proved  the  following: 

Theorem  3.10.  Let  X  =  {K,L),p  be  a  Bessel  system.  Let  X'  be  its  modified  adjoint.  Let  fl  be 
any  measurable  subset  of  IR*^,  such  that  {fl  -f  l}iec  pairwise  disjoint,  for  C  —  L,K.  Let  H  be 
the  subset  of  all  L20R'^)-functions  whose  Fourier  transform  is  supported  in  Q-hK.  Then,  there 
exist  unitary  maps 

V~  :  H  ^  £2(X'), 
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and 


Vl  :  ->  hiX), 


both  independent  of  ip,  such  that  the  identity 


(dmxr'/^n = 

holds  on  H.  Consequently,  given  any  Bessel  system  Y  =  {K,L),j,,  (denX)  ^TyT^  is  unitarily 
equivalent  on  H  to  the  restriction  ofTy/Tx'  ran  V~  C  hiX')- 

Remark.  Note  that  TJ^  maps  into  £2{X)  while  Ty  is  defined  on  e2(Y),  hence  TyT^  is,  formally, 
not  well-defined.  However,  the  two  index  sets  X,  Y  are  naturally  identified  with  the  set  {K,L), 

and  this  is  the  way  TyTx  should  be  interpreted. 

Theorem  3.10  falls  short  of  implying  directly  Theorem  2.2.  The  minor  reason  is  that  it  requires 
X  to  be  Bessel  (that  can  be  overcome  by  a  more  careful  analysis).  The  major  reason  is  that  we  do 
not  know  in  advance  that  the  various  norm  bounds  that  we  are  after  are  realized  on  spaces  H  of 
the  form  that  appear  in  the  theorem.  The  fact  that  the  relevant  bounds  are  realized,  indeed,  on 
such  spaces  H  is  an  immediate  consequence  of  Theorem  3.2,  which,  however,  leads  to  the  direct 
proof  of  Theorem  2.2,  presented  in  the  previous  subsection. 

Nevertheless,  Theorem  3.10  does  admit  some  applications.  For  example,  it  provides  the  missing 
proof  for  assertion  (b)  of  Theorem  2.2: 

Proof  of  (b)  of  Theorem  2.2.  If  X  is  not  fundamental,  then  Txf  =  0,  for  some  /  G  L2(IR‘^)\0. 
Let  X  =  and  let  be  any  subset  of  IR*^  such  that  (a):  /  does  not  vanish  identically 

on  A  :=  +  k,  (b):  The  K-shifts,  as  well  as  the  L-shifts,  of  are  pairwise  disjoint^  Without 

loss,  we  may  assume  /  to  be  supported  on  A  (otherwise,  we  define  g  &  L2  by  g  :=  /x,  with  x 
the  support  function  of  A,  and  apply  the  representation  of  in  terms  of  to  conclude  that 
T^g  =  0,  too).  Defining  H  in  Theorem  3.10  with  respect  to  the  present  fi,  the  theorem  provides 
US  with  the  relation 

Since  I4,  V~  are  partial  isometries,  this  readily  implies  that  Tx'  is  not  injective,  hence  so  is  Tx*. 
The  converse  is  obtained  by  a  similar  argument.  ^ 

A  stronger  version  of  Theorem  3.10  is  available  under  slightly  more  restrictive  conditions. 

Corollary  3.11.  Let  X  =  {K,L)^  be  a  Bessel  system,  X'  its  modified  adjoint.  Let  and  Ql 
be  fundamental  domains  for  the  lattices  K,  L,  respectively, 

(a)  If  the  L-shifts  of  are  pairwise  disjoint,  then  there  exist  unitary  maps  1^  :=  L2(1R  )  -> 
£2(X'),  and  V2,  such  that 

{denX)~^^^Tx  =  k^Tx'Vi. 

(b)  If  the  k -shifts  off}  L  are  pairwise  disjoint,  then  there  exist  unitary  maps  Vi  :  £2(X')  -t  ^2(1^'*), 
and  V2  such  that 

(denX)-^/^V2T^Vi=Tx>. 
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In  particular,  (a)  is  satisfied  by  all  univariate  systems  of  low-density,  and  (b)  is  satisfied  by  all 
univariate  systems  of  high-density. 

Proof.  For  the  proof  of  (a),  we  choose  :=  in  Theorem  3.10.  We  then  recognize  that 
~  S'lid  the  asserted  result  then  follows  from  that  theorem. 

’  For  (b),  we  choose  9,  :=  9l  in  Theorem  3.10.  The  construction  details  of  V~  that  precede  the 
theorem  then  imply  that  maps  ^  onto  £2{X').  We  then  take  Vi  :=  V~  and  V2  :=  V^. 

For  the  univariate  system  {2Tip7L,q7L)^,  we  may  choose  9^  :=  (0,1/p),  and  9l  =  (0, 9). 
Therefore,  one  of  the  intervals  is  included  in  the  other,  and  it  is  easy  to  relate  the  type  of  inclusion 
to  the  “right”  type  of  density.  ♦ 

Discussion.  The  last  corollary  applies  to  multivariate  systems,  such  as  sup-adjoint,  sub-adjoint  (cf. 
§4),  and  many  others,  but  certainly  there  are  multivariate  systems  that  satisfy  neither  condition. 
To  make  sure,  only  low-density  systems  can  satisfy  condition  (a),  and  only  high-density  can  be 
satisfy  condition  (b).  This  is  in  analogy  with  the  self-adjointness,  a  property  which  in  general 
implies,  and  in  the  univariate  is  also  implied  by,  the  denX  =  1  property. 

3.3.  The  Wexler-Raz  identity 

We  sidetrack  briefly  here  to  discuss  the  Wexler-Raz  identity  [WR],  whose  proof  in  [DL]  was  a 
key  ingredient  in  the  approach  there.  As  we  mention  in  the  introduction,  this  section  was  added 
to  the  present  article  only  after  we  became  aware  of  the  [DL]  work  and  its  details. 

The  intention  is  to  show  how  this  identity  is  related  to  our  fiberization  techniques,  or,  in  other 
words,  how  this  identity  reads  (hence  is  established)  in  terms  of  the  Gramian  matrices.  To  make 
sure  (as  should  be  the  case  with  any  identity,  and  in  contraist  with  norm  estimates),  we  do  not 
require  here  any  of  the  results  of  [RSI]  (such  as  Theorem  3.2),  and  hardly  any  of  the  results  here. 
Rather  though,  we  will  show  that  the  identity  reads  at  the  fiber  level  as 

vAu*  =  uA*v*, 

with  A  a  fiber  of  the  pre-Gramian  matrix,  and  u,v  some  £2-vectors  (of  the  right  “order”).  Such 
identity,  of  course,  holds,  but  only  after  one  verifies  that  all  sums  above  are  well-defined,  in  the 
sense  that  each  one  of  the  rows  and  columns  of  A  is  in  £2  (a  condition  which  is  self-evident  for  pre- 
Gramian  fiber  A),  that  Au*  as  well  as  vA  are  in  £2,  and  that  v{Au*)  =  {vA)u*.  For  the  satisfaction 
of  these  latter  requirements,  we  will  assume  that  all  our  WH  systems  are  Bessel,  and  will  invoke 
Theorem  3.2.3  of  [RSI]  which  rigorously  justifies  the  matrbc  manipulations  we  employ. 

Theorem  3.12  [WR],  [DL].  Let  X  =  {K,L)^,  Y  =  {K,L),i„  Z  =  {K,L)g  he  three  WH  systems. 
Assume  that  X*,  Y,  Z,  (and  Z*)  are  Bessel.  Then, 

(3.13)  TyT^<p  =  (denZ) 

Proof.  The  Bessel  assumption  on  the  systems  involved  ensures  us  that  both  sides  of  (3.13) 
are  well-defined  L2(Il‘^)-functions.  Let  r  (respectively,  r)  be  the  Fourier  transform  of  the  left  (resp. 
right)  hand  side  of  (3.13).  We  will  show  that  r  =  r,  a.e. 
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The  pre-Gramian  setup  detailed  at  the  beginning  of  §3.2  show  that,  for  a.e.  w  E  IR^, 

ri  ^ 

w+K 

with  V  :=  (j)\  Here,  we  wish  only  to  compute  r(u;),  hence  are  interested  only  at  the  {k  =  0)~row 

ly+if  ^ 

of  Jy{w).  That  row  consists  of  (cf.  the  third  display  of  §3.1)  u  :=  \K\  Altogether,  for 

a.e.  w, 

r{w)  - 

Analogous  computation  reveals  that,  for  a.e. 

r{w)  =  (denZ)|Zr^/Vj^*Hu  =  \K\^^^'^{denZ)^^^u^7^{w)v. 

Invoking  (3.1)  (with  X  there  replaced  by  Z*)^  we  conclude  that,  indeed,  r{w)  =  r{w).  ^ 

3.4.  Gramian  estimates 

Given  a  WH  system  A  =  (A:,L)$,  Theorem  3.2  provides  us  with  various  characterizations  of 
the  “basis  properties”  of  X,  Since  G{w)  and  G{w)  are  Hermitian,  we  may  obtain  useful  information 
by  studying  the  operator  spectrum  of  these  matrices.  Though  in  general  it  is  not  easy  to  compute 
the  spectral  radius  and  related  quantities  of  the  Gramian  and  dual  Gramian,  sufficient  conditions 
that  guarantee  the  boundedness  above  and  below  of  G{w)  can  be  easily  derived.  The  former  by 
various  standard  methods,  the  latter  by  “diagonal  dominance”  arguments.  Several  examples  of  this 
type  are  discussed  in  this  section.  The  arguments  that  lead  to  such  estimates  are  only  sketched.  In 
any  event,  these  are  straightforward  arguments  which  were  already  discussed  in  more  detail  (and 
in  a  more  general  context)  in  §1.6  of  [RSI]. 

For  the  simplicity  of  the  presentation,  we  will  consider  only  the  case  of  a  singleton  The 
extension  to  a  finite  $  is  mainly  notational. 

Results  analogous  to  some  of  the  statements  here  can  be  found  in  [Dl]  and  in  [T02];  more 
specific  pointers  are  given  in  the  sequel.  However,  it  is  worth  mentioning  a  difference  between  the 
results  in  the  aforementioned  references  and  the  ones  to  be  detailed  here.  That  difference  may 
be  easily  understood  by  the  following  illustration;  suppose  that  {fi)iei  ^  collection  of  functions 
defined  on  some  common  domain  A.  Then,  in  general,  the  condition 

Bi  :=  ^supl/il  <  00, 
ia 


is  stronger  than  the  condition 

B2  :=  sup^l/il  <  00. 
iei 

Furthermore,  one  usually  has  B2  <  Bi,  The  results  and  estimates  of  [Dl]  and  [T02]  are  based  on 
the  finiteness  and  “smallness”  of  quantities  defined  as  Bi  above,  while  our  results  below  are  based 
on  quantities  defined  similarly  to  the  number  B2  above,  and  thus  are  finer  than  their  univariate 
counterparts  in  the  literature. 
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Our  first  objective  is  the  derivation  of  upper  and  lower  bounds  for  ||r||  that  are  verifiable  by 
a  mere  inspection  of  the  entries  of  G  and  G.  For  that  purpose,  we  introduce  the  following  map, 
which  will  be  referred  to  as  the  ^transform  (M”  for  Ambiguity): 

A:=Ac-.  L2{JR‘^)  Li(lR‘^  Af{s,  •)  := 

lec 

Here,  C  is  some  fixed  lattice,  and  is  a  fundamental  domain  for  that  lattice. 

In  terms  of  the  .4-transform,  the  (f,  Z')-entry  of  the  Gramian  G{w)  of  the  WH  system  {K,  L)^ 
is 

The  (fc,  A:')-entry  of  the  dual  Gramian  G{'w)  is 

\K\-^AL(p{k-k',w  +  k'). 

These  observations  lead  us  to  the  following  estimates: 

Theorem  3.14.  Let  X  be  the  WH  system  {K,L)ip,  ip  G  L2(1R‘^).  Then: 

(a)  X  is  Bessel  system  if  the  function 

f-=Yl 

keK 

is  essentially  bounded  and  only  if  the  function 

9‘=i^ 

keK 

is  essentially  bounded,  and  we  have 

lialli,„(E*)  <  <  ll/lli,„(R-)- 

(b)  X  is  Bessel  system  if  the  function 

leL 

is  essentially  bounded  and  only  if  the  function 

l€L 

is  essentially  bounded,  and  we  have 

ll^illi^(ni-)  <  \K\\\Txf  <  ll/ilk^(iR-). 

Proof.  In  view  of  Theorem  3.2,  computing  HTx  ||  (and  thereby  verifying  the  Bessel  property 
of  X)  is  equivalent  to  computing  the  oo-norm  of  either  Q  of  Q*.  Since  G{w)  is  Hermitian,  we  may 
bound  its  2-norm  Q{w)  from  above  by  its  1-norm.  With  /i  as  in  (b),  ||/i||loo  easily  proved 

to  coincide  with  the  essential  supremum  (over  w)  of  the  1-norms  of  G{w).  A  similar  argument, 
with  G  replacing  G,  leads  to  the  estimate  that  involves  /. 

Again,  since  G{w)  is  non- negative,  we  can  also  bound  Q{w)  from  below  by  the  ^2"norm  of  each 
row  of  it.  This  leads  to  the  estimate  that  involves  gi .  Replacing  G  by  G  and  repeating  the  same 
idea,  leads  to  the  estimate  that  involves  g.  ^ 
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Discussion.  It  seems  instructive  to  pause  and  compare  the  different  estimates  provided  by  (a)  and 
(b)  in  the  above  theorem.  For  example,  the  upper  bound  of  (a)  is  majorised  by  the  sup-norm  of 


(3.15)  [^(- +k+ i)(p{- + 01 = iy(- + 01  y~i  i^(' + fc  -H  01, 

k&K,l&L  keK 

while  the  upper  bound  of  (b)  is  majorised  by 

(3.16)  |y(-  +  k)\ |^(-  +  k  +  l)\. 

keK 

Assume,  for  example,  that  L  C  K.  Then,  (3.15)  becomes 

(3.17)  |^(•  +  fc)l  + 

keK 

As  to  (3.16),  identifying  K /L  with  some  finite  A  C  K,  that  expression  can  be  written  as 

(3.18)  +  '  + 

6eA  l£L 

While  pointwise  the  two  estimates  are  not  comparable,  it  is  easily  seen  that  the  sup-norm  of  (3.18) 
is  <  the  sup-norm  of  (3.17).  Since  the  second  estimate  is  based  on  the  Gramian,  and  the  example 
just  studied  is  of  a  low-density  system,  this  suggests  that  the  Gramian  upper  bound  estimates  ((b) 
in  the  theorem)  are  better  (i.e.,  tighter)  for  low-density  systems,  hence  that  the  dual  Gramian 
upper  bound  estimates  (i.e.,  (a)  of  the  theorem)  yield  better  results  for  high-density  systems. 

In  any  event,  the  above  analysis  shows  that  T  is  bounded  whenever  (p  decays  fast  enough,  or 
equivalently,  whenever  p  is  smooth  enough: 

Corollary  3.19.  Let  X  be  a  PWH  system  generated  by  p.  Then  X  is  a  Bessel  system  if 

a€C 

is  essentially  boundedj  for  £  =  L,  K. 

Proof.  From  the  above  discussion,  it  follows  that 

keK  keK 

Theorem  3.14  (a)  then  implies  the  claim.  ^ 
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Discussion  cont’ed.  Analogous  estimates  are  available  in  terms  of  ip  (rather  than  (p)  with  the 
lattices  L,  K  replaced  by  their  duals  K,  L.  Consequently,  A  is  a  Bessel  system  if  either  ip  is 
sufficiently  smooth  or  decays  sufficiently  fast. 

Literature  Discussion.  We  compare  the  /-estimate  in  (a)  of  Theorem  3.14  to  Theorem  2.6  of 
[Dl]  and  Theorem  2  of  [T02].  Theorem  2.6  contains  an  estimate  of  the  upper  frame  bound  for 
a  univariate  X  in  terms  of  ip  rather  than  ip.  After  making  that  necessary  switch,  we  find  that  it 
employs  quantities  of  the  form  ||  YIi^l  ||oo-  However,  these  norms  are  not  summed  up 

directly,  but  rather  cire  grouped  into  pairs,  with  respect  to  which  a  geometric  average  is  computed, 
hence  altogether  Theorem  2.6  is  not  comparable  to  our  results  here.  As  to  Theorem  2  of  [T02], 
rewriting  it  in  terms  of  ip  rather  than  p,  its  sufficient  condition  for  the  Bessel  property  of  the 
univariate  X  is  the  finiteness  of  Xligz,  with  ai  the  Fourier  coefficients  of  That 

boundedness  condition  clearly  implies  the  boundedness  of  /  in  Theorem  3.14. 

The  derivation  of  sufficient  conditions  for  A  to  be  a  Bessel  system  is  much  simpler  than 
guaranteeing  A  to  be  a  frame.  Unless  we  are  willing  to  settle  for  conditions  for  a  fundamental 
frame  A  or  a  Riesz  basis  A,  it  is  virtually  impossible  to  derive  feasible  sufficient  conditions  for 
frames  that  are  based  on  the  magnitude  of  various  entries  of  G  or  G.  However,  in  order  for  A 
to  be  a  fundamental  frame,  G{w)  should  be  invertible,  and  that  can  be  observed  in  case  G(w)  is 
diagonally  dominant.  Sufficient  conditions  of  this  type  are  discussed  in  the  next  theorem. 

Theorem  3.20.  Let  X  —  {K,  L)^  be  a  Bessel  PWH  system. 

(a)  If  the  function 

f  :=  AiipiO,  •)  -  \ALP{k,  •)! 

fceK\o 

is  positive  a.e.  and  is  (essentially)  bounded  away  from  0,  then  X  is  a  fundamental  frame  for 
L2(IR‘*)  and 

(3.21)  \\T*-Y  <  I^11|1//IIl,.(ir<')- 

(b)  If  the  function 

g  ■.=  )  -  E 

l^L\0 

is  positive  a.e.,  and  is  (essentially)  bounded  away  from  0,  then  X  is  a  Riesz  basis,  and 

||T-if  <  |A|||l/5|U^(n^«). 

Proof.  Part  (b)  follows  from  (a)  and  the  duality  principle:  to  prove  that  A  is  Riesz  basis, 
it  suffices,  by  that  principle,  to  prove  that  A*  is  a  fundamental  frame.  A  conversion  of  the  condition 
in  (a)  from  (A,  to  (L,K)^  then  yields  the  condition  in  (b). 

We  now  prove  (a):  by  Theorem  3.2,  we  need  to  show  that  each  Q{w)*~  is  bounded.  Note  that 
the  entries  of  each  row  of  the  dual  Gramian  G{w)  are  of  the  form 

\K\~^ALip{k,w'),  keK, 
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for  an  appropriately  chosen  E  IR^.  The  condition 

Al${0,w)-  ^  \AL^{k,w)\  >  e,  a.e.  «;  e 

keK\0 

thus  implies  that,  for  almost  every  w  G  IR^  and  for  every  row  of  G{w),  the  diagonal  entiy 
AL(piO,w')  is  £  greater  than  the  £i-normj3f  the  other  entries  of  this  row.  Thus,  as  soon  as  this 
condition  holds,  the  self-adjoint  operator  G{w)  is  diagonally  dominant,  hence  invertible,  for  almost 
every  w,  and  g*-{w)  =  llG(u>)-i||  <  By  Theorem  3.2,  X  is  a  fundamental  frame,  and  the 

estimate  (3.21)  holds.  ^ 

Remark.  A  slightly  different  variant  of  the  above  theorem  is  as  follows:  if  /  in  the  theorem  is  >  e 
on  the  support  o/.4l<?(0,  •),  then  X  is  a  frame  (not  necessarily  fundamental)  and  the  estimate 

(3.21)  holds  with  ||l//||Loo(n)  replacing  ||1//11z,oo(ir‘')- 

Literature  Discussion  cont’ed.  Theorem  2.5  of  [Dl]  provides  a  sufficient  condition  for  a  uni¬ 
variate  system  to  be  a  fundamental  frame.  The  condition  is  based  on  the  positivity  of  (a  slightly 
coarser  expression  than)  the  quantity 

inf.4.L^(0,  •)  “  ^  ||.4L^(fc)  Olloo) 

keK\0 

and  thus  the  condition  given  in  Theorem  3.20  is  somewhat  better.  In  any  event,  both  conditions 
are  based  on  very  coarse  estimates.  The  comparison  of  Theorem  3.20  to  the  lower  bound  estimate 
provided  in  Theorem  7  of  [T02]  is  similar. 

As  stated  several  times  before,  the  “basis”  (K,  L)^  satisfies  the  same  properties  as  its  Fourier 
transform  set  {L,K)^.  This  means  that  all  the  above  estimates  are  valid  with  (p  replacing  (p  and 
AT,  L  interchanging  roles. 

The  above  Theorem  3.14  and  Theorem  3.20  are  “pointwise”  in  the  sense  that  we  hold  the 
lattices  K  and  L  fixed.  In  the  last  part  of  our  Gramian  analysis,  we  refine  certain  “asymptotic 
results”  from  [Dl]  (see  also  [D2]).  The  idea  is  as  follows:  suppose  that  we  fix  L  and  p,  but  vary 
the  volume  of  K.  As  K  becomes  denser,  X  is  more  likely  to  become  a  fi:ame.  This  prediction  is 
nicely  reflected  in  the  dual  Gramian:  as  K  becomes  denser,  K  become  sparser;  thus,  while  the 
diagonal  entry  of  the  (fc=0)-row  of  \K\G{w)  (we  are  allowed  to  consider  only  such  row;  cf.  the  proof 
of  Theorem  3.20)  i.e.  Alp{^i  w),  is  independent  of  the  choice  of  K,  the  1-norm  of  the  off-diagonal 
entries  ALp{k,w)  is  sensitive  to  the  sparsity  of  K.  Specifically,  if  p  decays  reasonably  fast  (i.e., 
if  p  is  smooth  enough),  we  expect  Alp{s,u>)  to  decay  to  0  in  a  controllable  way  as  |s|  oo. 
Analogously,  if  p  is  nicely  localized,  we  may  switch  from  {K,  L)^  to  (L,  K)^  (hence  will  oversample 
the  modulations  rather  than  the  shifts).  So,  roughly  speaking,  under  some  basic  conditions  on 
the  smoothness  (decay)  of  p,  we  should  be  able  to  obtain  a  frame  from  the  Bessel  system  X  by 
increasing  the  density  of  K  {L,  respectively).  A  special  result  in  this  direction  follows. 
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Theorem  3.22.  Let  X  :=  {K,  L)^  be  a  Bessel  system.  DeGne 

^{s)  :=  ^l(s)  :=  ML^(sr)llLoo(iR‘')- 


Assume  that 

(a)  The  sequence  summable  (i.e.,  lies  in  Ci{K)). 

(b)  The  shifts  {tp{-  -  form  a  Riesz  basis  (respectively,  a  frame)  with  Riesz  (frame)  bounds 

A  and  B. 

Then,  for  all  sufGciently  large  integer  n,  the  PWH  system  X„  :=  {K/n,L),p  is  a  fundamental 
frame  (frame,  respectively).  Furthermore, 

lim„_,oo^^^^  =  denXj5. 
n 

lim„^oo  =  denX  A, 

n 

Proof.  We  prove  only  the  case  when  F  :=  {v?(-  -  0}/^^  is  a  Riesz  basis.  It  is  well- 
known  (cf.  e.g.  [RSI],  [BW]),  that  F  being  a  Riesz  basis  (or  frame)  is  equivalent  here  to  the 
following  two  conditions:  (i):  /?(0)  <  oo  (which  is  equivalent  to  F  being  a  Bessel  system),  and 
(ii):  S(0)  :=  ess  inf  ^1,^(0,  •)  >  0.  (for  frame:  S(0)  :=  essinfo.(^)  •)  >  0,  where  o{F)  is 

the  spectrum  of  F-,  cf.  [RSI]  for  details).  It  is  furthermore  known  that  the  frame  bounds  of  F  are 
B  =  WT^r  =  m/\L\  and  A  =JTf-X^  =  3(0)/|Z|. 

Define  Kn  :=  K/n.  Then,  Kn  =  nK,  and  hence,  since  is  summable,  c„  :=  ->  0 

as  n-^  oo.  By  Theorem  3.14, 


jji;^  <  =  |L||fr:r'(||TFf  +  jZr'cn)  ->  (denX)IITFf  =  (denX)B. 

n  n|A„| 

Further,  the  same  theorem  provides  the  simpler  estimate 


^i^^>-^  =  (denX)R, 
n  nlRTnl 


and  the  first  claim  is  thus  established.  ' 

We  now  prove  the  other  claim.  Choosing  n  large  enough,  we  can  guarantee  that  (S(0)— c„)  >  0. 
Theorem  3.20  then  implies  that  is  a  frame  and  that 

--i* . -  >  ->  iZlIA'r^A  =  (denX)A. 

n  n\Kn\ 

Again,  the  converse  inequality  is  straightforward.  By  the  definition  of  S(0)  we  can  find,  for  any 
e  >  0,  a  set  of  positive  measure  fig  such  that  Ai,ip{0,w)  <  3(0)  +  e  for  every  w  £  fie.  Since 
\Kn\~^AL<p{^,  w)  is  a  diagonal  element  of  Gx„(io),  (regardless  of  the  value  of  n),  we  conclude  that 
^*Xn  —  d*  on  By  letting  e  — 0,  we  obtain 

n||Tx„rMP  =  >  l^l/S(0)  =  ((denX)A)-^ 
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The  theorem  suggests  the  construction  of  “snug  frames”  (Daubechies  terminology,  meaning 
that  the  frame  is  “almost”  tight)  by  “going  to  the  limit”  with  a  smooth  orthonormal  system;  one 
starts  with  an  orthonormal  system  F  :=  {(^(-  —  V  decays  fast  enougl^,  one  expects 

the  function  p2L‘‘  ^e  summable.  This  means  that  for  almost  every  lattice  K,  is 

summable.  By  taking  X  :=  {Kln,L)^  for  sufficiently  large  n,  one  obtains  an  “almost  tight”  frame 
for  L2(1R-‘^)-  If  F  is  “merely”  tight  frame,  the  same  procedure  applies,  only  that  the  fundamentality 

of  X  is  lost. 

As  alluded  to  before,  an  alternative  construction  applies  to  functions  which  decay  nicely,  but, 
say,  are  not  very  smooth.  In  that  case,  one  requires  {^(-  —  Q!)}a6z;‘'  orthonormal,  defines  the 
function  /?  with  respect  to  cp  instead  of  0,  and  proceeds  by  oversampling  the  modulations  rather 
than  the  translations.  Again,  the  construction  in  [Dl]  can  be  viewed  as  a  particular  instance  of 
that  strategy.  Indeed,  in  case  d  =  1,  the  following  result  is  intimately  related  to  Theorems  2.5  and 
2.6  of  [Dl]. 

Theorem  3.23.  Let  X  :=  {K,L)^  be  a  Bessel  system.  Define 

0{s)  :=  Pxis)  :=  ||.4i(:<^(s,-)||Loo- 


Assume  that  _ 

(a)  The  sequence  Is  summable  (Le.,  lies  in  h{L)). 

(b)  The  K-modulations  of  (p  form  a  Riesz  basis  (a  frame),  with  Riesz  (frame)  bounds  A  and  B. 
Then,  for  all  sufficiently  large  integer  n,  the  PWH  system  :=  (K,L/n)^  is  a  frame  (which  is 
fundamental  in  the  Riesz  case).  Furthermore, 

limn_oo^^^^  =  denAB, 

n 


-111-2 


lim- 


rx.r'ii 


n— )>oo 


n 


=  denXA, 


We  also  remark  that  some  improvements  are  available  in  case  (respectively,  /3x)  satires 
slightly  better  decay  conditions  than  the  mere  assumption  ^  ^i(^)  (respectively,  Pk  ^ 

For  example,  if  /?l  (respectively,  Pk)  is  majorized  by  a  radially  symmetric  non-increasing  summable 
function  /?'  (respectively,  )9'),  then  it  is  easy  to  see  that  for  any  lattice  iiT,  (respectively,  L)  and  for 
sufficiently  large  positive  (not  necessarily  integer)  n,  the  system  (iiT/n,  L)^p  (respectively,  (iiT,  L/n)(^) 
is  a  frame. 
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4.  Zak  transform  analysis 


4.1.  Multivariate  Zak  transform 

The  Gramian  estimates  that  were  derived  in  the  previous  section  are  crude  in  the  sense  that, 
generally  speaking,  they  fail  to  take  into  account  possible  cancellations  in  the  application  of  the 
operators  G{w)  and  G{w).  However,  there  are  cases  when  the  Gramian  G(w)  and  the  dual  Gramian 
G(w)  are  intimately  connected  with  a  finite  collection  of  convolution  operators,  and,  in  such  cases, 
required  quantities  like  0(w),  Q(w)~,  etc.,  can  be  more  accurately  computed  via  the  symbols  of 
the  underlying  convolution  operators.  This  section  is  devoted  to  the  development  of  this  approach. 
The  symbols  of  the  relevant  convolution  operators  are,  in  turn,  related  to  the  Zak  transform  of  the 
generating  function  (p,  and  this  is  the  reason  for  the  title  chosen  for  the  present  section.  As  already 
mentioned  in  the  introduction,  the  reference  [ZZ]  contains,  in  the  univariate  setup,  results  which 
are  equivalent  to  some  of  the  results  listed  in  this  section. 

We  define  the  multivariate  analog  of  the  Zak  transform  as  follows.  First,  let  K  be  some 
lattice  in  IR*^.  Let  I  be  any  linear  isomorphism  between  K  and  its  dual  K  (e.g.,  choose  a  basis  for 
K,  map  it  to  its  dual  basis  in  K  and  extend  the  map  by  linearity).  Then 

Z^f{w,t)  :=  f{w  +  k)eikit),  f  e  L20R^). 

k^K 

When  considering  f  as  defined  on  9.k  x  Ox,  Z^  is  an  isometry  from  a  dense  subspace  of 
L2{JR'^)  into  a  dense  subspace  of  I-2(Ox  x  Ox)  (with  the  latter  being  normalized  appropriately), 
hence  extends  to  an  isometry  between  these  two  spaces.  This  fact  is  well-known  in  the  univariate 
case  (cf.  p.  976  of  [Dl])  and  extends  to  the  multivariate  case  with  no  difficulty.  The  details  of  the 
map  I  will  be  insignificant,  and  we  will  suppress  it  and  write 

Z^f{w,t)  =  Y  + 

keK 

In  an  earlier  draft  of  the  present  paper,  the  course  of  analysis  of  this  section  was  based  on  the 
duality  principle,  but  later  we  found  a  way  to  present  the  results  without  invoking  that  principle; 
still  the  notions  that  were  introduced  in  the  context  of  the  duality  principle  are  very  helpful  here  as 
well.  For  example,  the  simplest  case  in  the  present  analysis  is  the  case  of  a  self-adjoint  X,  since  then 
each  J{w),  G{w),  and  G{w)  is  simply  a  convolution  operator.  In  fact,  the  same  is  true  for  G{w) 
even  for  a  sup-adjoint  X,  and  even  when  the  singleton  (p  is  replaced  by  finitely  many  generators. 
In  this  case,  the  dual  Gramian  operator  can  be  compressed  into  a  single  function,  a  function  which, 
in  a  univariate  setup,  was  already  put  into  good  use  in  [Dl]  (being  the  way  Daubechies  derives 
her  “exact  bounds”,  cf.  p.  982  of  [Dl]).  We  therefore  start  our  discussion  with  an  analysis  of  a 
sup-adjoint  X  =  {K,L)^. 

4.2.  Analysis  of  sup-adjoint  and  sub-adjoint  WH  systems 

By  definition,  if  {K,  L)  is  sup-adjoint,  then  K  C  L.  We  set 

r 
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for  the  finite  group  L/K,  The  same  notation,  F,  will  be  used  in  the  sequel  to  denote  any  set  of 
representers  (in  L)  of  this  group.  Introducing  the  set 

r4>:={MV: 


we  see  that 

(4.1)  {K,L)^=^{K,K)r,^ 

This  means  that  every  sup-adjoint  system  X  can  be  viewed  as  a  WH  system  (iiT,  L)^,  whose 
corresponding  lattice  pair  {K^L)  is  self-adjoint  Therefore,  without  loss,  we  may  assume  to  this 
end  that  X  =  (iT,  _ 

We  now  study  the  pre-Gramian  of  such  X,  Here,  the  rows  of  Jx  are  indexed  by  $  x  iF. 
We  organize  in  row-blocks  ((/?,  iF),  G  and  denote  the  corresponding  block  by  J*.  Then,  we 
observe  that,  for  any  fixed  w  G  IR*^,  the  rows  of  J^{w)  comprise  the  set  of  all  iF-shifts  of  sequence 

a^,«,  :K  ^C:kh^  \K\-^^‘^<p{w  +  k). 

Consequently  (with  possible  some  minor  re- indexing),  the  map 

'•  ^2{K)  -)■  i2{K)  :  c  1-^  J*{w)c 


is  the  convolution 


This  implies  that  the  entire  map  Jxi'^)  ^  vector- valued  convolution  operator: 


(4.2)  •  ^2{K)  ^2(^  ^  K)  •  ^  ^  (^V?,iy  * 

It  easily  follows  from  that  that  the  operator  Jx  (w;)  can  be  identified  with  the  convolution  operator 
Jx{w)  :  ^2{^  xK)  ^  e2{K)  :  h4  ^  a^,wi—)  *  Cp- 


Altogether,  we  obtain  the  following. 

Proposition  4.3.  Let  X  be  the  self-adjoint  WH  system  {K,K)^,  and  $  C  L2(1R‘^)  be  finite. 
Then,  for  each  to  e  IR*^  the  dual  Gramian  operator  G{w)  is  a  convolution  operator  of  the  form 

G{w)  :  f2{K)  ^2{^)  •  ^  ^  ^  ^  *  ®v>,tu(  ■)  * 

Theorem  3.2  reduces  the  study  of  various  properties  of  X  to  a  corresponding  study  of  the  fiber 
operators  {G{w))^,  In  the  present  case,  G{w)  is  a  convolution  operator,  and  therefore  all  properties 
and  quantities  in  question  may  be  studied  via  the  symbol  of  the  convolution,  i.e.,  its  Fourier  series. 
In  the  present  case,  we  need  to  compute  the  symbol  of 


As  one  easily  computes,  the  symbol  in  question  is  the  function 


1-^1  ^  S  I  XI  + 

k£K 

This  symbol  can  be  written  in  terms  of  the  Zak  transform  of  ^  as 

Before  summarizing,  we  convert  the  above  expression  from  the  self-adjoint  case  to  the  sup-adjoint 
case.  In  the  latter  situation,  the  set  $  is  replaced  by  the  set  r$,  T  :=  LfK  (cf.  (4.1))  and  therefore 
the  symbol  becomes  _ 

'y^L/ K 

Our  observation  concerning  sup-adjoint  systems  X  can  be  now  summarized  as  follows. 

Corollary  4.4.  Let  X  be  a  sup-adjoint  system  {K,  L)$.  Then,  for  a.e.  w  E  IR'^,  the  dual  Gramian 
operator  is  a  convolution  operator  with  symbol 

Zx{w,-):=\K\-^Y1  X  1^  +  7,  OP- 

j^L/K 

Standard  properties  of  convolution  operators  then  imply  the  following: 

Corollary  4.5.  Let  X  =  {K,  L)$  be  a  sup-adjoint  WH  system.  Then,  in  the  notations  of  Theorem 
3.2,  and  with  Zx  defined  as  above,  we  have 

(a)  Q*x{w)  =  \\Zx{w,-)\\l^. 

(b)  G*x{w)  =  \\llZj:{w,-)\\ll. 

(c)  G*x~\{w)  =  l|l/-2x(t«,-)lli^(,^„),  with  the  support  of  Zxiw,-). 

Combining  the  last  theorem  with  Theorem  3.2,  we  obtain  the  following  result  concerning  sup- 
adjoint  X.  The  univariate  counterpart  of  (a)  and  (b)  of  this  result  is  known  and  can  be  found  e.g., 
in  [Dl]. 

Theorem  4.6.  Let  X  =  iK,L)i,  $  C  La  hnite,  K  C  L.  Set  T  :=  L/K,  and  denote 

zx:= 


Then  the  following  is  true: 

(a)  X  is  a  Bessel  system  if  and  only  if  Zx  G  Loo(lR‘^  x  IR‘^)-  Furthermore,  l|Tx|P  =  li^xlP  = 
W^xWloo- 

(b)  Assume  X  is  Bessel.  Then  X  is  a  fundamental  frame  if  and  only  if  Zx  is  bounded  below  (away 

from  zero).  Furthermore,  =  ||l/-2xllz,„o- 

(c)  Assume  that  X  is  Bessel.  Then  X  is  a  frame  if  and  only  if  Zx  is  bounded  below  on  its  support. 

Furthermore,  ||Txi"Mr^  =  with  a  the  relevant  support. 
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The  above  analysis  also  implies  that,  for  a  sup-adjoint  X  =  the  operator  is 

realized  on  the  “Zak  transform  domain”  as  the  multiplication 

(4.7)  L2  9  /  •->  f' 

To  make  sure,  deriving  the  last  theorem  directly  from  that  observation  is  quite  straightforward, 
in  stark  contrast  with  the  amount  of  work  required  for  proving  Theorem  3.2  (which  we  invoked 
in  our  approach  here).  However,  having  already  Theorem  3.2  at  our  disposition,  our  argument 
contributes  further  to  the  understanding  of  the  connection  between  the  Gramian  analysis  and  the 
Zah  transform  analysis. 

Since,  in  our  Gramian  analysis,  we  did  not  state  results  concerning  the  fundamentality  of  X, 
we  had  no  result  of  this  kind  to  transport  into  the  Zak  transform  analysis.  Nevertheless,  it  is  easy  to 
observe  from  (4.7)  that  X  is  fundamental  if  and  only  if  Zx  vanishes  almost  nowhere;  further,  the  L2- 
functions  whose  Zak  transform  is  supported  on  the  complement  of  supp  Zx  comprise  the  orthogonal 
complement  of  X.  Also,  since  the  frame  operator  TxT^  is  realized  on  the  Zak  transform  domain 
as  multiplication  by  Zx,  the  dual  frame  operator  is,  necessarily,  division  by  the  same  function, 
implying  that  the  Zak  transform  of  the  generating  functions  of  the  dual  frame  is  given  by 

Z^ip/Zx, 

a  well-known  phenomenon  for  univariate  fundamental  frames.  Here,  0/0  :=  0.  More  discussion 
along  these  lines  is  given  near  the  end  of  this  section,  when  we  compare  sub-adjoint  systems  to  PSI 
spaces  (in  this  regard,  compare,  also,  (b)  in  the  corollary  below  to  Theorem  2.2.16  of  [RSI]). 

Corollary  4.8.  Let  X  be  the  sup-adjoint  WH system  (K,  L)^,  with  dual  Zak  function  Zx ,  whose 

support  is  cr{X).  Then:  „ 

(a)  The  closed  span  of  X  consists  of  all  L2-functions  whose  Zak  transform  Z  is  supported  in 

- 

(b)  Assuming  X  to  be  a  frame,  its  dual  frame  is  of  the  form  (K,  L)ip,  where,  for  each  'll)  Z  'ip 
is  supported  on  cr{X)  and  defined  there  by 

Z^i)  :=  {Z^ip)/Zx- 

In  case  $  is  a  singleton,  the  duality  principle  (i.e.  Theorem  2.2)  allows  us  to  obtain  several 
useful  results  with  respect  to  sub-adjoint  systems,  i.e.,  systems  of  the  form  {K,L)(p  where  L  G  K 
(or,  equivalently,  PWH  systems  X  that  satisfy  X  C  X*).  Here,  X*  is  sup-adjoint,  and  the  function 
Zx*<,  written  directly  in  terms  of  X  (rather  than  in  terms  of  X*)  is 

seA 

with 

A  :=  KjL. 

The  duality  principle,  when  combined  with  Theorem  4.6,  then  connects  the  basic  properties  of  the 
sub-adjoint  X  to  the  function 

Zx  :=  (denXjZx.  =  \K\-^  X) 

(56A 

The  implied  result  (which  can  certainly  be  derived  directly,  too)  is  then  the  following: 
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Theorem  4.9.  Let  X  be  a  PWH  sub-adjoint  system  (K,L)^.  DeHne 

S6A.=K/L 

Then  the  following  is  true. 

(a)  X  is  a  Bessel  system  if  and  only  if  Zx  €  Loo  Furthermore,  ||Tx|P  =  =  I1-2a^I1loo- 

(b)  Assume  X  is  Bessel.  Then  X  is  a  Riesz  basis  if  and  only  if  Zx  is  essentially  bounded  below 

away  from  zero.  Furthermore,  =  ||1/'2^x||l«,- 

(c)  Assume  X  is  Bessel.  Then  X  is  a  frame  if  and  only  if  Zx  is  essentially  bounded  below  on  its 
support  a.  Furthermore,  ||rx|“^ll"^  =  \\l-l ZxWlooW)- 

Discussion:  Sub-adjoint  systems  as  the  WH  analogue  of  PSI  spaces.  A  PSI  (=principal 
shift-invariant)  space  S{<j))  is,  by  definition,  the  La-closure  of  the  shifts  of  the  La-function  (f).  The 
study  of  F  :=  ^  ^  potential  basis  for  S{4>)  is  the  simplest  among  all  studies  of  shift- 

invariant  bases  (cf.  [RSI],  [BW]  and  the  references  therein).  Up  to  some  technicalities,  a  function 
/  G  S{4>)  is  identified  by  the  restriction  of  its  Fourier  transform  to  [0  . .  2tt]^,  hence  the  rest  of  its 
Fourier  transform  values  are  somewhat  “redundant”  information.  Closely  related  to  that  is  the  fact 
that  each  of  the  “fiber  spaces”  (obtained  by  restricting  the  Fourier  transform  of  S{<f))  to  a  lattice  of 
the  form  «; -I- 27rK‘^)  is  either  one-  or  zero-dimensional.  The  Gramian  matrix  here  is  of  order  1,  and 
its  single  entry  is  the  function  [^,^]  :=  Y,a&2n7L<‘  \^i'  +  which  simply  computes  the  square  of 
the  norm  of  all  fibers  of  <^.  It  follows  that,  on  the  frequency  domain,  the  relevant  operator  TxTx 
is  realized  as  the  multiplication  operator 

^2(K*^)  9  c  1-4  [^, (^]c  G  L2(IR*^). 

The  support  of  [$,  $]  is  the  spectrum  of  S{(l>)  (i.e.,  it  is  the  set  of  all  one-dimensional  fibers)  and 
is  useful  for  finding  ker  Tx  and  its  orthogonal  complement. 

Analogously,  in  the  case  of  a  sub-adjoint  WH  X  :=  {K,L),p,  a  function  /  in  the  closed  span 
of  X  is  essentially  identified  by  the  restriction  of  its  Fourier  transform  to  fl~  -t-  L.  The  fibers  here 
are  indexed  by  x  and  each  fiber  space  is  an  one-  or  zero-dimensional  subspace  of  ^2(^) 
(spanned  by  the  evaluation  of  the  pre-Zak  vector  at  the  corresponding  iw,t)).  The  operator  TxTx 
is  realized,  on  “the  Zak  transform  domain”  as  multiplication  by  Zx-  We  define  the  WH  spectrum 
of  X  to  be  the  support  a{X)  C  x  of  Zx-  The  generator  V*  of  the  dual  frame  of  X  can  be 
found  exactly  as  in  the  PSI  case  (cf.  (2.2.13)  in  [RSI]),  that  is,  it  can  be  computed  fiber  by  fiber: 

:=  Z^!plZx, 

with  the  division  carried  out  only  on  the  WH  spectrum  of  X  (more  precisely,  on  the  L-periodization 
of  that  spectrum;  one  should  be  somewhat  careful  here,  since  Z^(p  has  as  its  fundamental 

domain,  while  Zx  has  fl~y.fl.  as  its  fundamental  domain.  Thus,  these  expressions  should  be 
extended  to  an  L-periodic  expression  in  the  first  argument  before  the  division  can  take  place).  The 
analogy  goes  on  and  on.  For  example,  a  suitable  unitary  transformation  which  maps  onto 

L2{fl^  X  n^)  maps  kerTx  onto  L2{flj^\a{X),flL)  (compare  with  (2.2.8)  and  its  previous  display 
in  [RSI]). 


31 


Our  final  discussion  in  this  section  concerns  with  oversampling  and  undersampling  of  WH 
systems:  Each  sup-adjoint  system  can  be  realized  as  the  result  of  oversampling  a  self-adjoint  one, 
and  each  sub-adjoint  system  can  be  realized  as  the  result  of  undersampling  a  self-adjoint  system. 
Let  X  =  {K,K)^  be  the  self-adjoint  system,  Y  =  {K,L)^  its  sup-adjoint  counterpart  (that  is 
K  C  L)^  and  Y*  the  corresponding  sub-adjoint.  If  X  is  a,  fundamental  frame,  then  it  is  more  or  less 
trivial  to  see  that  Y  is  so,  too,  and  hence  (by  the  duality  principle,  say)  that  y  is  a  Riesz  basis. 
It  is  much  less  obvious  to  see  that  if  X  is  a  frame  (not  necessary  fundamental),  then  both  Y  and 
y*  are  frames,  too;  it  follows,  however,  from  the  analysis  of  the  present  subsection.  In  contrast, 
the  converse  does  not  follow:  thus,  while  X  can  be  viewed  as  a  multi-sampled  version  of  Y  ,  the 
possible  frame  property  of  y *  may  be  lost  during  that  process.  In  summary,  such  oversampling 
does  not  preserve  in  general  the  frame  property, 

4.3.  Zak  transform  analysis  of  compressible  WH  systems 

In  the  previous  subsection,  it  was  shown  that,  for  a  sup-adjoint  and  sub-adjoint  WH  systems 
X,  the  study  of  (sup- adjoint  case)  or  Tx  (sub-adjoint  case)  is  reduced  to  studying  the  behaviour 
of  a  single  function  {Zx  for  sup-adjoint,  Zx  for  sub-adjoint).  Though  the  same  does  not  hold  for 
other  WH  systems,  there  are  more  general  situations  when  the  analysis  of  X  can  be  reduced  to 
the  study  of  finite- order  Hermitian  matrices.  In  fact,  the  situation  here  is  analogous  to  that  that 
occurs  in  the  study  of  shift- invariant  spaces.  While,  as  discussed  in  the  previous  subsection,  the 
PSI  space  and  the  sub- adjoint  WH  space  have  one-dimensional  fiber  spaces,  the  more  general  FSI 
(=finitely  generated  shift-invariant)  spaces  are  similar  to  the  present  compressible  WH  systems: 
the  fiberization  of  both  lead  to  fiber  spaces  of  finite  dimension.  In  the  latter  case,  each  fiber  space 

Is  spanned  by  the  columns  of  PZx{nj^  f),  hence  is  considered  a  subspace  of  ^2(^)  (s^^  below  for 
definitions  and  details).  The  critical  information  required  is  that  basic  operations  (such  as  finding 
image  and  kernel  of  operators,  describing  orthogonal  projectors,  computing  dual  system,  etc.)  can 
be  performed  fiber  by  fiber.  For  FSI  spaces,  this  was  done  in  §2  of  [RSI];  the  techniques,  however, 
apply  to  general  fiberization  with  finite-dimensional  fibers,  hence  to  our  present  case  of  interest. 

We  introduce  a  new  type  of  “inner  product” ,  the  Zak  product^  of  which  the  function  Zx  that 
was  used  in  the  context  of  the  sup-adjoint  case  is  a  special  “diagonal”  case.  We  will  then  form 
finite-order  matrices  (the  Zak  matrix  and  the  dual  Zak  matrix)  whose  entries  are  such  Zak  products. 
These  matrices  may  replace  the  infinite  order  Gramian  matrices,  and  we  reach  that  reduction  by 
employing  either  of  the  two  arguments:  (a)  direct  adaptation  of  the  FSI  techniques  of  [RSI]  to 
WH  setup  here;  or  (b)  adaptation  of  the  FSI  techniques  in  order  to  fiberize  further  the  (already 
partially  fiberized)  Gramian  operator,  and  then  invoking  Theorem  3.2.  Our  remark  from  the  last 
section  concerning  the  latter  approach  is  valid  here  as  well:  it  is  much  easier  to  derive  the  results 
here  directly  than  proving  Theorem  3.2  (as  one  observes  by  comparing  the  analysis  of  FSI  spaces 
and  general  SI  spaces  in  [RSI]). 

We  call  a  lattice  pair  (X,  L)  compressible  if  the  group  X  fl  L  has  a  finite  index  in  L  (equiva¬ 
lently,  in  X),  which  is  the  case  exactly  when  X  fl  L  is  d-dimensional.  We  define  then  two  quotient 
groups: 

r:-L/(Lnx), 
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and 


A:=  kl{Lr\k). 

We  also  think  about  T,  A  as  any  set  of  representers  for  the  above  groups.  We  refer  to  the  order  of  T 
as  the  compression  factor  of  X,  and  to  the  order  of  A  is  the  decomposition  factor  of  X.  (The 
terminology  reflects  the  fact  that  we  are  interested  primarily  in  high-density  systems  X.  For  low- 
density  systems,  the  notions  of  compression  and  decompositions  factors  should  be  interchanged). 
Note  that  X  is  self-adjoint  if  and  only  if  both  F  and  A  are  trivial,  X  is  sup-adjoint  if  and  only  if  A 
is  trivial,  and  X  is  sub-adjoint  if  and  only  if  F  is  trivial.  Note  also  that  the  notations  F ,  A  retain 
their  meaning  from  the  last  section  in  case  of  a  sup/sub-adjoint  system. 

Example  .  Suppose  that  X  =  {2Trp7Z,q7L)^.  Then,  knL  =  7Llpn  qTL.  Clearly,  k  (Mis  non¬ 
trivial  iff  pq  is  rational.  Since  pq  =  (denJf)  we  conclude  that  a  univuriats  is  compressible 

iff  its  density  is  rational. 

In  the  case  of  a  sup-adjoint  X,  we  identified  the  pre-Gramian  with  a  vector- valued  convolution 
operator.  It  is  not  hard  to  see  (as  we  are  going  to  show),  that  in  the  present  case  the  pre-Gramian 
is  a  matrix-valued  convolution  operator. 

Here  are  the  details:  we  let  X  =  {K,L)^,  ^  finite,  be  a  compressible  WH  system.  The  rows 
of  the  pre-Gramian  are  indexed  by  L  x  $  =  (X  fl  L)  x  F  x  $,  and  we  re-organize  them  in 
blocks  {{K  n  L),j,  (/?),  ^  €  $,  7  G  F.  The  columns  of  Jx  are  indexed  K  =  {K  H  L)  x  A,  which  we 
also  organize  in  blocks  {{k  n  L),  S),  6  e  A.  This  induces  a  block  re-organization  of  the  entire 
with  a  typical  block  indexed  by  ((7,  v’),^)-  As  in  the  sup-adjoint  case,  one  observes  that  the  block 
*^(7  v>)  <5  ^  convolution  operator 

o-y,(p,s,w  •  n  L)  — t  t2{K  n  L) :  c  i->  *  c. 

Here,  the  sequence  is  defined  by 

(4.10)  a7,v),«,ti;(j)  :=  +  d  +  j). 

The  operator  Jx(ta)  acts  from  =  f-2{{K  fl  L)  x  A)  into  i2{L  x  $)  =  £2{{K  fl  L)  x  F  x  $), 

and,  by  the  above,  acts,  indeed,  as  a  matrix-valued  convolution  operator: 

Jxiw)  :  (ci)ieA  ^  (XI  ®7,v>.i.«'  *  C5)7€r,v’e$- 

«GA 

This  observation  allows  us  to  replace  the  pre-Gramian  by  the  matrix  of  the  symbols  of  the  various 
0'y,(p,5,w  As  in  the  sup-adjoint  case,  the  symbol  of  the  sequence  is  the  Zak  transform 

Thus,  we  obtained  a  representation  of  the  pre-Gramian  Jxi^)  as  a  finite  order  matrix  whose  rows 
are  indexed  by  F  x  $,  whose  columns  are  indexed  by  A,  and  whose  structure  is  as  follows: 

PZ*x{w,-)  :=  +  7  + 
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It  is  quite  obvious  that  the  analogous  representer  of  Jx(^)  is  the  adjoint  •)  of  the  above 

PZ^(iu,  •)•  We  refer  to  both  of  these  matrices  as  pre-Zak  matrices.  The  spectrum  a{X)  of  X 
is  the  set  of  all  {w^t)  where  PZxi'io^t)  is  not  the  zero  matrix. 

The  development  now  becomes  quite  transparent.  Using  the  pre-Zak  matrices,  we  construct 
finite-order  analogs  of  the  Gramian  and  the  dual  Gramian.  The  analog  of  G{w)  =  JxJxi'^) 
square  non-negative  (F  x  $)-order  matrix  *),  whose  ((7, (^),  (7',^))-entry  is 


where,  for  any  lattice  £,  finite  set  A  C  IR^,  and  /,p  €  1^2  ?  Z£(/,5)  is  the  Zak  product  of  /  and 
g  defined  as 

(4.11)  Za(/,5)  :=  Zi{f,g)  :=  + 

SeA 

In  analogy  with  the  Gramian  matrix,  we  call  the  matrix  2x  the  Zak  matrix.  The  analog  of 
the  dual  Gramian  is  the  dual  Zak  matrix  obtained  as  Zx  •—  PZxPZ*x-  Thus,  it  is  a  square 
non-negative  A-order  matrix,  whose  (5, 5')-entry  is 

\K\-^  if). 

Example.  Let  m,n  be  positive  integers,  g.c.d.(m,n)  =  1.  Let  X  be  the  univariate  WH  system 
{2Tr7L/m,n7L)^  (with  $  finite,  otherwise  arbitrary).  Then  K  =  mTL,  hence  K  O  L  =  mnTL.  We 
may  thus  take  F  =  L/{K  f]  L)  —  (n,  2n,  ...,mn),  and  A  =  K/{K  fl  L)  =  (m, 27n,  ...,nm).  The 
relevant  Zak  transform  here  is  Z  :=  which  can  be  defined  as 

Zf{w,t)  :=  ^  fiw  +  mnj)e2T,j/mn{t)- 
j€7L 


The  pre-Zak  matrix  PZx  has  the  form 

(m/27r)^/^  {Z0{-  +  jn  +  j'm, 

(where  rows  are  indexed  by  j'  and  columns  by  Consequently,  the  dual  Zak  matrix  is  an 

n  X  n  matrix  whose  (j,f)-entry  is 

771 

(m/27r)  ^  ^  (Z(p)(w  ■+  in  +  jm,t){Z(fi){w  +  in  + 


The  Zak  matrix,  say  in  case  $  is  a  singleton  (/?,  is  an  m  x  m  matrix  whose  (j,  j')-entry  is 

n 

(m/27r)  Y^{Z(p){w  -b  im  +  jn,t){Z(p){w  -H  zm  -b  j'n,t). 
i=l 
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Finally,  the  following  variant  of  the  duality  principle,  though  is  not  an  issue  here,  arises  in  a  natm-al 
way:  interchanging  between  m  and  n  results  at  a  similar  interchange  between  the  Zak  matrix  and 
the  dual  Zak  matrix.  At  the  same  time,  such  an  interchange  amounts  to  passing  from  X  to  its 
adjoint  X*.  Indeed,  the  fact  that  passing  to  the  adjoint  amounts,  up  to  a  multiplicative  constant, 
to  an  interchange  between  the  two  basic  matrices  is  the  essence  of  the  duality  principle. 

Remark.  Zibulski  and  Zeevi,  [ZZ],  as  well  as,  implicitly,  Daubechies,  [Dl],  employ  the  Zak 
transform  in  order  to  decompose  univariate  systems  of  rational  density.  In  particular,  [ZZ]  employs 
two  matrices,  S  and  G,  which  correspond  to,  yet  look  quite  different  from,  our  dual  Zak  matrix,  and 
pre-Zak  matrix.  However,  the  difference  can  be  simply  attributed  to  a  different  choice  of  a  basis  for 
4(A)  or  4(r).  Our  Zak  matrices  here  are  related  to  the  standard  bases  for  these  spaces  (made  of 
functions  of  one-point-support).  On  the  other  hand,  the  [ZZj-matrices  are  related  to  choosing  the 
characters  of  the  dual  group  as  the  relevant  basis.  Since  both  choices  are  orthonormal,  our  dual 
Zak  matrix  here  is  unitarily  equivalent  to  the  S  matrix  of  [ZZj. 

In  order  to  convert  the  analysis  of  the  “basis  set”  X  from  the  infinite-order  Gramian  matrices  to 
the  finite  order  Zak  matrices,  we  need  to  know  how  to  relate  the  functions  that  appear  in  Theorem 
3.2  {Q,  G*,  etc.)  to  the  Zak  matrices.  However,  as  was  already  explained  before,  the  study  of  the 
connection  between  G{'w)  and  Zxiw,-)  here  is  entirely  analogous  to  the  connection,  in  the  case 
of  an  FSI  set  X,  between  the  standard  infinite-order  Gramian  {{x,y))x,y&x  and  the  compressed 
Gramians  (G(u;)),„giRd.  These  latter  relations  were  analysed  in  detail  in  several  references,  with 
the  most  comprehensive  results  contained  in  [BDR]  and  [RSlj.  In  particular,  the  arguments  used 
to  establish  Theorem  2.3.6  of  [RSI],  when  transported  to  the  new  setup  here,  lead  to  the  following 
result. 

Lemma  4.12.  Given  a  compressible  WH  system  X  :=  (K,  L)$  and  its  associated  Zak  matrix  Zx, 
consider,  for  fixed  w,t  E  IR*^,  the  matrix  Zxiw,  t)  as  a  map  from  4(r  x  $)  into  itself.  Let  Cx{w,  i) 
be  the  norm  of  this  map,  ^x{w,t)  the  norm  of  the  inverse  (defined  as  oo  if  Zx{w,t)  is  singular), 
and  (^i(w,  t)  the  norm  of  the  pseudo-inverse.  Then,  in  the  notations  of  Theorem  3.2,  the  following 
relations  hold: 

(a)  For  a.e.  w  6  IR**, 

Gxi'w)  =  ||Cx(«^,*)I|loo- 

(b)  For  a.e.  w  G  IR*^, 

^x(^)  ~  IICx(^>  ■)  11^00  • 

(c)  For  a.e.  w  6  IR*^, 

with  the  support  of  Cx(u^>  •)• 

Invoking,  thus.  Theorem  3.2,  we  obtain  the  following  characterization  of  the  “basis”  properties 
of  the  compressible  X  in  terms  of  the  Zak  matrix. 

Theorem  4.13.  Let  X  =  (K,L)^  be  compressible,  and  let  Zx  be  its  associated  Zak  matrix. 
Consider,  for  each  {w,t)  E  IR*^  x  IR*^,  the  matrix  Zx(‘w,t)  as  a  map  from  4(r  x  $)  into  itself,  and 
let  Cx{u!,t),  (^(w,t)  and  (^i(w,t)  be  the  associated  norm  functions  defined  in  Lemma  4.12.  Then: 
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(a)  X  is  a  Bessel  system  if  and  only  if  Cx  ^  ^  IR*^).  Furthermore,  llTxH  11^x11 

IlCxIUoo-  -r^-  r  n  ff, 

(b)  Assume  X  is  a  Bessel  system.  Then  X  is  a  Riesz  basis  if  and  only  ifQx  ^  I'urttiermore, 

K'f  =  IICxIIloo-  _  .. 

(c)  Assume  X  is  Bessel.  Then  X  is  a  frame  if  and  only  if  Cx\  ^  Loo(criX)),  witha(X)  the  spectrum 

ofX.  Furthermore,  ||Txi“^lP  =  IICx|llioo(<r(x))- 


While  the  Zak  matrix  is  useful  for  computing  the  important  functions  Q,  Q  and  associated 
with  the  Gramian  of  X,  the  dual  Zak  matrix  is  usefuHor  computing  the  complementary  quantities 
Q*,  Q*~  and  Q*~  associated  with  the  dual  Gramian  G.  For  that,  we  consider,  for  each  {w,t)  €  IR  , 

the  matrix  Zx{w,t)  as  a  map  from  ^2(A)  into  itself,  and  denote  by  Cx(^>*)  (“^x  Cx]  ) 

the  norm-function  (respectively,  the  inverse-norm  function,  and  the  pseudo-inverse-norm  function) 
of  Zx.  Again,  the  arguments  used  in  Theorem  2.3.6  of  [RSI]  apply  here,  and  lead  to  a  result 
analogous  to  Lemma  4.12,  that  connects  now  the  C*-functions  to  the  ^*-functions.  This,  in  view 
of  Theorem  3.2,  leads  to  following  theorem. 


Theorem  4.14.  Let  X  be  a  compressible  WH  system  (K,  L)$.  Let  Zx  be  its  associated  dual  Zak 
matrix,  and  let  Cx,  Cx"  and  Cx~\  be  the  corresponding  norm- functions.  Then: 

(a)  X  is  a  Bessel  system  if  and  only  if  Cx  ^  Loo(lR^  x  IR  ).  Furthermore,  ||Tx||  —  11^x11 


(b) 

(c) 


Assume  X  is  Bessel  Then  X  is  a  fundamental  frame  if  and  only  if  Cx  ^ 

\\Trr  =  \\Cx-\\L^-  .  r  . 

Assume  X  is  Bessel  Then  X  is  a  frame  if  and  only  if*  Cx|  ^ 

ll^x|  ^11^  =  IICx|  IUoo(o-(-^))* 


FurthermorOj 

Furthermore, 


Remark.  For  a  univariate  X,  parts  (a)  and  (b)  of  the  above  result,  as  well  as  part  (b)  of  the 
following  corollary,  are  due  to  [ZZ]  (Theorem  2  and  Proposition  1  there).  Though  Theorem  2  in 
[ZZ]  is  stated  without  proof,  the  supporting  discussion  there  provides  an  almost  complete  argument. 
The  only  missing  ingredient  in  the  approach  of  [ZZ]  is  its  sought-for  unitary  diagonalization  (14), 
which,  incidentally,  is  proved  in  [RSI]  as  Lemma  2.3.5. 

Corollary  4.15.  Let  X  be  a  compressible  WH  system  with  Zak  matrix  Zx  and  dual  Zak  matrix 

Zx  •  Then:  _ 

(a)  X  is  a  tight  frame  if  and  only  if,  after  normalizing  the  generator  of  X,  either  Zx  or  Zx  (or 
both)  is  an  orthogonal  projector  a.e.  (equivalently,  if  the  matrix  spectrum  of  almost  every 
Eber  of  the  Zak  or  dual  Zak  matrix  consists  of  at  most  0  and  1). 

(b)  X  is  a  fundamental  tight  frame  if  and  only  if  Zx  is  a  constant  multiple  of  the  identity  a.e.  on 

IR*^  X  IR*^. 

(c)  X  is  orthonormal  if  and  only  if  Zx  is  the  identity  a.e.  on  IR  x  IR  . 

Proof.  Since  A  is  a  tight  frame  ilf  ||Tx|l  =  ||Zx^^ll’"^  ^e  obtain  from  Theorem  4.13  that 
the  tightness  of  X  is  equivalent  to  the  equality  HCxIUoo  ~  IICx|IUooi  which  clearly  amounts  to  the 
fact  that  Zx{w,t)  has  only  one  eigenvalue,  independent  of  {w,t),  other  than  0.  Tins  eigenvalue 
is  the  normalization  constant  appearing  in  (a).  The  statement  with  respect  to  Zx  is  proved 
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similarly,  by  invoking  Theorem  4.14,  and  recalling  that  tightness  is  also  equivalent  to  the  equality 

IKII  =  Kf'll-'- 

Statement  (b)  follows  from  (a),  when  combined  with  (b)  of  Theorem  4.14. 

Statement  (c)  follows  easily  from  the  last  part  of  Theorem  3.2.  A 

Theorem  4.14  ,  in  particular,  shows  that,  in  order  for  X  to  be  a  fundamental  frame,  almost 
all  dual  Zak  matrices  Zx{w,t)  must  be  invertible.  In  fact,  as  is  already  established  in  Theorem 
1  of  [ZZ],  the  frame  property  of  X  is  irrelevant  here:  By  invoking  Lemma  2.3.5  of  [RSI],  we 
may  unitarily  diagonalize  Zx,  using  a  unitary  matrix  with  measurable  entries.  Thus,  if  the  fibers 
{Zx{w,  <)}(«,, t)  of  Zx  are  singular  on  a  set  of  positive  measure,  one  can  easily  construct  a  function 

/  such  that  T^T^f  =  0  (the  fibers  of  /,  4-  S,t))seA  should  each  lie  in  the  kernel  of 

Zx(w,  t)  for  a.e.  {w,  t),  and  the  only  reason  we  need  Lemma  2.3.5  here  is  in  order  to  synthesize  the 
fiber  kernels  into  a  measurable  function).  However,  Zx{w,t)  is  non-singular  iff  the  pre-Gramian 
PZx{w,t)  is  of  rank  #A,  the  later  being  possible  only  if 

#A  <  #r  X  ##. 

Since  ^A  is  the  index  of  if  D  L  in  X,  and  #r  is  the  index  of  if  fl  L  in  L,  we  see  that 

#A  ^  \KnL\  JL|  ^  M  ^  1 

#r  \K\  |ifnL|  \k\  den(if,L)' 

Defining,  naturally,  the  density  of  X  as 

X  den  (if,  L), 

we  obtain  the  following  result.  The  univariate  case  of  this  result  can  be  found  in  [Dl],  and  was 
proved  there  by  a  similar  argument. 

Corollary  4.16.  Let  X  =  (if,  L)$  be  a  compressible  WH  system.  Then  X  is  fundamental  only  if 
denX  >  1. 

Finally,  we  consider  the  problem  of  computing  the  dual  frame  of  X  =  {K,L)qf.  We  first 
study  this  problem  when  X  is  a  Riesz  basis,  say,  for  H  C  L2(1R‘^).  Let  R  be  a  map  from  $  to 
L2(IR‘^).  In  order  for  Y  :=  (if,  L)r$  to  be  bi-orthogonal  with  X,  we  need  TyTx  to  be  the  identity, 
hence  that  PZyPZx  is  the  identity  matrix  a.e.  Furthermore,  for  a  dual  basis  Y,  the  column- 
space  of  (almost)  each  PZy(it;,t)  should  be  the  same  as  that  of  PZx{w,f)  (follows  from  the  fact 
that  kerTj^  =  kerT^).  Thus,  finding  PZY{w,t)  here  is  a  standard  finite-dimensional  least  squares 
problem,  and  we  obtain  that 

(4.17)  PZY{w,t)  =  Zx^{w,t)PZx{w,t),  for  a.e.  (w;,f); 

(compare  with  Theorem  2.4.7  of  [RSI]). 
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The  situation  is  even  simpler  for  a  fundamental  frame:  with  the  understanding^hat  inverting 
Zx  amounts  to  the  pointwise  inversion  of  each  of  its  fibers  Zx{w,  t),  we  realize  that  represents 
the  inverse  of  the  dual  frame  operator  TxTx-  Thus  the  Zak  t^nsform  representation  of  the 
generators  R$  of  the  frame  dual  to  X  can  be  found  by  applying  Z^^  to  the  representation  of 
Since  the  action  TxTx  on  /  is  represented  in  the  form 

Zx  (•+^,  •))<56A) 

we  see  that  (cf.  (29)  of  [ZZ]) 

(4.18)  PZY{w,t)  =  Zx'^{w,i)PZx{w,t),  a.e.. 

In  the  case  of  a  non-fundamental  frame  X,  the  same  argument  can  be  employed,  but,  alas,  the 
pseudo- inverse  of  Zx  should  be  computed.  In  the  special  case  when  Zx{w,t)  is  either  non-singular 
or  0,  the  computation  of  the  pseudo-inverse  is  avoided,  and  this  explains  why  we  do  not  have  to 
tackle  this  problem  in  the  sup-adjoint  case. 

4.4.  Zak  transform  estimates 

The  Zak  transform  analysis  had  led  us  to  the  decomposition  of  the  operators  TxTx  '^x'^xi 
for  a  compressible  WH  set,  into  the  fibers  {Zx{w^t))^  t^^d  and  {Zxiw,t))^  i;^^d.  In  §3,  we  had 
exploited  the  Gramian  fiberization  in  order  to  estimate  the  frame  (Riesz)  bounds.  In  the  same 
manner,  we  may  exploit  here  the  Zak  transform  fiberization  is  order  to  derive  analogous  estimates. 
These  estimates  are  collected  below  without  further  explanation  (other  then  pointing  to  their  §3 
counterparts).  One  should  observe  that  despite  of  the  seemingly  close  relation  between  the  Gramian 
estimates  and  the  Zak  transform  estimates,  there  is  practically  substantial  difference  between  the 
two:  the  information  required  for  the  Gramian  estimates  is  readily  available  (i.e.,  values  of  either  tp 
or  ip\  we  are  tacitly  assuming  that  p  is  explicitly  known,  where  “explicitly”  might  mean,  e.g.,  that 
p  ox  (p  are  given  analytically);  at  the  same  time,  the  Gramian  estimates  are  crude.  In  contrast, 
the  Zak  transform  estimates  require  finer  information  (certain  Zak  transforms  and  subsequently 
Zak  products),  but  award  us  with  better  estimates  (particularly  when  the  relevant  Zak  matrix  is 
of  small  order). 

As  we  did  in  §3.2,  we  assume  here  that  X  is  the  principal  compressible  (K,  L)^.  The  compress¬ 
ibility  of  X  is  essential  (otherwise,  the  Zak  matrices  remain  of  infinite  order).  The  extension  to 
non-principal  systems  is  primarily  notational.  All  Zak  transforms  of  this  subsection  are  computed 
with  respect  to  the  lattice  K  r\  L.  In  fact,  the  only  Zak  transform  which  is  required  is  that  of  p, 

hence  we  set  _ 

g  :=  Z^^^p. 

As  before,  T  and  A  are  the  quotient  groups 

r  =  L/{KnL),  A  =  Kl{KnL). 

It  is,  perhaps,  worth  noting  that  the  Zak  matrix  and  the  dual  Zak  matrix  are  invariant  under 
shifts  by  A:  G  XnL  (in  both  variables);  this  observation  entitles  us  to  represent  any  7  G  T  or  ^  G  A 
in  any  convenient  way,  as  well  as  to  switch  between  representers,  as  the  arguments  for  proving  the 
estimates  below  require. 


38 


Theorem  4.19.  Let  X  =  {K,L)^  be  a  compressible  PWH  system.  Let  g,  T,  and  A  be  as  in 
the  second  paragraph  of  this  subsection.  Then  X  is  a  Bessel  system  if  and  only  if  g  is  essentially 
bounded.  Furthermore,  let  5a, 7  be  the  function 

5a,7  :=  ^  5('  +  7  +  ')9{'  +  ')) 

seA 

and  let  gr,s  be  the  function 

gr,s  +  7, -M-  +  7, •)• 

7€r 


Then: 

(a) 

=  l|Tif  <|Ar|-‘llEl9A,,lllL. 

7er 

7€r 

(b) 

I^T'IKE 

5eA 

=  l|rif<|Jir|-‘|iEl9r,j|||L„ 

SeA 

Proof.  We  note  that  \K\~^l'^g{w,  t)  is  one  of  the  entries  of  the  pre-Zak  matrix  PZx  (w,  t), 
and  that,  further,  every  entry  of  PZx(w,t)  is,  up  to  a  unit  multiplicative  constant,  of  the  form 
IKl~^^^g(w',t')  for  (possibly  other)  Since  the  pre-Zak  matrix  is  of  finite  order,  this  proves, 

in  view  of  Theorem  4.13,  the  fact  that  X  is  Bessel  iff  g  is  bounded. 

The  bounds  asserted  in  parts  (a)  and  (b)  are  proved  by  an  argument  analogous  to  that  used 
in  the  proof  of  Theorem  3.14,  after  observing  that  7)767  comprise  the  entries  of  a  typical 

row  of  Zx,  while  |iirl“^(5r,«)5eA  comprise  the  entries  of  a  typical  row  of  Zx-  ♦ 

“Diagonal  dominance”  arguments  similar  to  the  ones  employed  in  the  proof  of  Theorem  3.20, 
when  combined  with  the  observations  just  made  in  the  proof  of  Theorem  4.19,  lead  to  the  following 
sufficient  conditions  for  X  being  a  Riesz  basis  or  a  fundamental  frame. 

Theorem  4.20.  Notations  as  in  Theorem  4.19.  Let  X  be  a  compressible  Bessel  system;  then  the 
following  hold. 

(a)  X  is  Riesz  basis  if  the  equality 

9a,o  ~  1pa,7|  ^  £ 

7er\o 

holds  for  some  positive  e.  Furthermore,  if  this  is  the  case,  we  also  have 

||Tx-‘f  <  \K\le. 

(b)  X  is  a  fundamental  frame  if  the  equality 

9r,o  -  IPr.il  >  £ 

seA\o 

holds  for  some  positive  e.  Furthermore,  if  this  is  the  case,  we  also  have 
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